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Abstract 



Recently-generated long strong-coupling series for the two-point Green's func- 
tions of asymptotically free O(A^) lattice a models are analyzed, focusing 
on the evaluation of dimensionless renormalization-group invariant ratios of 
physical quantities and applying resummation techniques to series in the in- 
verse temperature (3 and in the energy E. Square, triangular, and honeycomb 
lattices are considered, as a test of universality and in order to estimate sys- 
tematic errors. Large-A^ solutions are carefully studied in order to establish 
benchmarks for series coefficients and resummations. Scaling and universality 
are verified. All invariant ratios related to the large-distance properties of the 
two-point functions vary monotonically with N , departing from their large- A'^ 
values only by a few per mille even down to A'^ = 3. 

PACS numbers: 11.15.Me, 75.10 Hk, 11.10 Kk, 11.15 Pg. 
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I. INTRODUCTION 



The properties of physical systems in the vicinity of a critical point, such as critical 
exponents and amplitude ratios, can be extracted by a variety of methods, ranging from 
exact solutions to Monte Carlo simulations. 

In the absence of exact results, one of the most successful approaches is based on the 
investigation of the strong-coupling series expansion, which enjoys the property of a finite 
radius of convergence, often (but not necessarily) coinciding with the extent of the high- 
temperature phase. More generally, when no singular points occur on the real axis of the 
complex coupling plane, it is possible to exploit strong-coupling results even beyond the 
convergence radius by analytic continuations, which are based on appropriate resummation 
methods. Extending the length of the strong-coupling series and improving the accuracy of 
the resummations are therefore the two most compelling tasks within this approach to the 
study of the behavior of systems in the critical region. 

As part of an extended program of strong-coupling calculations we have recently com- 
puted an extended series expansion of all nontrivial two-point Green's functions 



for the nearest-neighbor lattice formulation of two-dimensional O(A^) a models on the square, 
triangular, and honeycomb lattices, respectively up to 21st, 15th, and 30th order in the 
strong-coupling expansion parameter (3. A complete presentation of our strong-coupling 
computations for 0{N) a models in two and three dimensions will appear in a forthcoming 
paper. A preliminary report of our calculations can be found in Ref. [|1|. 

The relevance of a better understanding of 2-d 0{N) a models cannot be overestimated. 
They appear in condensed matter literature as prototype models for critical phenomena 
that are essentially restricted to two-dimensional layers, including some instances of high-Tc 
superconductivity. Moreover, they can be employed as model field theories sharing some of 
the most peculiar features of four- dimensional gauge theories, such as asymptotic freedom 
and spontaneous mass generation. This last statement must however be qualified, since 
the above-mentioned properties, according to common lore, are possessed only by those 2-d 
0{N) models such that > 2. 

We focus here on these asymptotically free models, analyzing their strong-coupling ex- 
pansion in order to extract information that may be relevant to the description of their 
continuum limit (/5 oo), assuming /3c = oo to be the only singularity on the real axis. 
This hypothesis is favored by all numerical evidence as well as by the successful application 
of the extrapolation techniques that we shall discuss in the present paper. The analysis of 
our strong-coupling series for models with A^ > 2 is presented in Ref. 0. 

It is obviously quite hard to imagine that strong-coupling techniques may be really 
accurate in describing the divergent behavior of such quantities as the correlation length 
and the magnetic susceptibility. Nevertheless, as our calculations will explicitly confirm, the 
strong-coupling analysis may provide quite accurate continuum-limit estimates when applied 
directly to dimensionless, renormalization-group invariant ratios of physical quantities. Two 
basic ideas will make this statement more convincing. 

(i) For any dimensionless, renormalization-group invariant ratio R{P), when (3 is suffi- 
ciently large we may expect a behavior 




(1) 
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where R* is the fixed point (continuum) value and is the (diverging) correlation length. 
Hence a reasonable estimate of R* may be obtained at the values of f3 corresponding to 
large but finite correlation lengths, where the function R{f3) fiattens. This is essentially 
the same idea underlying Monte Carlo studies of asymptotically free theories, based on the 
identification of the so-called scaling region. 

(ii) On physical grounds, it is understandable that /? is not necessarily the most conve- 
nient variable to parameterize phenomena occuring around (3 = oo. An interesting alterna- 
tive is based on the observation that the strong-coupling series of the internal energy 

E = /3 + 0(/33) (3) 

may be inverted to give /9 as a series in E. This series may be substituted into other strong- 
coupling expansions, obtaining expressions for physical quantities as power series in E. It 
might now be easier to reach the continuum limit, since it now occurs at a finite value of 
the expansion variable, i.e., E ^ 1. 

We hope to convince the reader that, by exploiting these ideas, state-of-the-art strong- 
coupling calculations can be made at least as accurate as the best Monte Carlo simulations 
presently available, when applied to dimensionless renormalization-group invariant quanti- 
ties. 

We must stress that the analysis of the strong-coupling series calculated on different 
lattices offers a possibility of testing universality, and, on the other side, once universality 
is assumed, it represents a further check for possible systematic errors and allows their 
quantitative estimate; this estimate is usually a difficult task in strong-coupling extrapolation 
methods such as those based on Pade approximants and their generalizations. 

Our physical intuition of the behavior of 0(iV) models is strongly guided by our knowl- 
edge of their large- behavior, and by the evidence of a very weak dependence on N of the 
dimensionless ratios. In order to extend our understanding to those lattices that have not 
till now received a systematic treatment, and also in order to establish a benchmark for the 
strong-coupling analysis, we decided to start our presentation with a detailed comparative 
study of the large- limit of various lattices, in the nearest-neighbor formulation. To the 
best of our knowledge, only the large- A^ solution on the square lattice was already known 
explicitly [^. 

The paper is organized as follows: 

In Sec. |T|we present the large- A^ limit solution of O(A^) a models on the square, triangu- 
lar and honeycomb lattices, in the nearest-neighbor formulation, calculating several physical 
quantities and showing explicitly the expected universality properties. The triangular- and 
honeycomb-lattice results are original, and possess some intrinsic reasons of interest. How- 



ever, readers willing to focus on square-lattice results are advised to jump to Sec. |T| after 
reading Subs. [11 A| and [11 B|, where the notation is fixed. 



Sec. m is devoted to a detailed analysis of the available strong-coupling series of G{x) 
and other physical quantities on the square, triangular, and honeycomb lattices. Most of 
the results we shall show there concern the A^ = 3 model. The basic motivation for this 
choice lies in the observation that all dependence in A^ is monotonic between 3 and oo; hence 
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the discussion of higher- results would be only a boring repetition of the considerations 
presented here. The reader not interested in the analysis of triangular and honeycomb 
lattices may skip most of the discussion, by focusing on Subs. [Ill B| , where further definitions 
are introduced and the square-lattice series are analyzed, and on Subs. [Ill E| , where all 
conclusions are drawn. 

Apps. and |B| provide the derivation and the technical details of the large- calculations 
on the triangular and honeycomb lattices respectively. We present as well the calculation of 
the A-parameters. 

App. y is a study of the complex temperature singularities of the N = oo partition 
functions on the triangular and honeycomb lattices. 

In Apps. 0, and |^ we present, for selected values of A^, the strong-coupling series of 
some relevant quantities on the square, triangular, and honeycomb lattice respectively. 



II. THE LARGE- AT LIMIT OF LATTICE 0{N) a MODELS 

A. The large- AT saddle point equation 

The nearest-neighbor lattice formulations on square, triangular and honeycomb lattices 
are defined by the action 

Sl = -Np J2 ^^i ■ ^^r, s^-S:^ = 1, (4) 

links 

where s* is a A^-component vector, the sum is performed over all links of the lattice and 
xi^Xr indicate the sites at the ends of each link. The coordination number is c = 4, 6,3 
respectively for the square, triangular and honeycomb lattice. The lattice spacing a, which 
represents the length unit, is defined to be the length of a link. The volume per site is then 
Vs = 1, \/3/2, 3-\/3/4 (in unit of a?) respectively for the square, triangular, and honeycomb 
lattice. 

Straightforward calculations show that the correct continuum limit of O(A^) a models, 

5* = — y (fxd^s{x) ■ d^s{x), s{x) ■ s{x) = 1, (5) 
is obtained by identifying 

respectively for the square, triangular and honeycomb lattice. Notice that 

\ = t/3 = ^ (7) 

c 

is the distance between nearest-neighbor sites of the dual lattice in unit of the lattice spacing 
a. 
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When the number of field components N per site goes to infinity, one can use a saddle 
point equation to evaluate the partition function. Replacing the constraint s'J' = 1 by a 
Fourier integral over a conjugate variable a^, we write the partition function as 



Zoc Yl dsx^isl - 1) exp N(3 ^ 4, ■ 4^ 

X links 



links 



(8) 



Integrating out the (p variables we arrive at the expression 



Z oc y dax exp — iox — Tr In i? 



(9) 



where 



and /S.xy is a generalized Laplacian operator, such that 

^ i'Pxi - (Pxr f = <Px^xy(i)y (H) 

links x,y 

The large- iV limit solution is obtained from the variational equation with respect to a^- 
Looking for a translation invariant solution we set 

iax = Y^- (12) 



The matrix R then becomes 



Rxy — ^ [ ^xy ~l" ZVg^xy] , (l^) 

and the saddle point equation is written as 

1 = lim —TyR-\ (14) 

where Ng is the number of sites. 

The large- iV fundamental two-point Green's function is obtained by 

G{x-y)=R-l. (15) 

In order to calculate the trace of -R~^, the easiest procedure consists in Fourier trans- 
forming the operator R. Such transformation is straightforward on lattices, like square and 
triangular lattices, whose sites are related by a translation group, and in these cases it yields 
the diagonalization of the matrix Rxy. The honeycomb lattice, not possessing a full trans- 
lation symmetry, presents some complications. In this case a partial diagonalization of Rxy 
can be achieved following the procedure outlined in Ref. H]. 
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B. The square lattice 



Turning to the momentum space the variational equation becomes 

^ _ r, _ r d'^k 1 _ 1 
t~^~ J-n (27r)2 k^ + z ~ 2^ 



7 = = r T^rr^ = ^P^m iPsiz)) , (16) 



where 



1 



-1 



p,{z) = {l + -z) , (17) 



and K is the complete integral of the first kind. 
Let's define the moments of G{x) 



m,,^Y.i^'yG{x). (18) 



Straightforward calculations lead to the following results 



t , . 

X = mo = -, (19) 



z 



il - ^ ^ = j, (20) 



m| 1 z 



Xm^ 4 V 16/ 

Notice that in the large- A'^ limit the renormalization constant of the fundamental field is 
Z = t. M is a renormalization-group invariant quantity. 

The mass-gap should be extracted from the long distance behavior of the two-point 
Green's function, which is also related to the imaginary momentum singularity of the Fourier 
transform of G{x). In the absence of a strict rotation invariance, one actually may define 
different estimators of the mass-gap having the same continuum limit. On the square lattice 
one may consider /ig and /i^ obtained respectively by the equations 

G~^{pi = iids,P2 = 0) = 0, 

lis and /Xd determine respectively the long distance behavior of the side and diagonal wall- wall 
correlations constructed with G{x). In generahzed Gaussian models, such as the large- A?" 
limit of 0{N) models, it turns out convenient to define the following quantities 

= 4 ( cosh^ - 1 ) . (23) 
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In the continuum limit 

Ms Md 



1, (24) 

/^s /id 

therefore Ms and Md may be also used as estimators of the mass-gap. 
In the large- iV limit 

= Mi = z = Ml. (25) 

The rotational invariance of G{x) at large distance, (i ^ ^, is checked by the ratios fis/f^d- 
Using the above results one can evaluate the scaling violation terms: 

/^d y21n(^v^+/rT|i) 48 23040 ^ ^ ^ ^ 

Another test of scaling is provided by the ratio 



^In (^ + ,[^]=l-^z + ^z^ + 0(z^). 
V 2 V 4; 24 640 ^ ^ 



Mg 

The internal energy can be easily calculated obtaining 



(27) 



E = (4 ■ = R-'.+, = 1 - ^ + i- (28) 



Therefore 



^E((^V.-4)^) = ^ - |, (29) 

where the term proportional to z is related to the condensate T of the trace of the energy- 
momentum tensor 

^d,s{x)-d,s{x). (30) 

In the large- iV limit 

m = (31) 



therefore from the expression (^8]) we deduce 

T _ 1 
M| ~ 4^' 



(32) 



Another interesting quantity which can be evaluated in the large-A^ limit is the zero- 
momentum four-point renormalized coupling constant, defined by 



7 



where 



X4= (^0 ■ 4 4 ■ ^z)c- (34) 

x,y,z 

is zero in the large- hmit, where the theory is Gaussian-hke and thus X4 = 0. Its value 
in the continuum limit 

can be also evaluated in the large- expansion of the continuum formulation of the O(A^) 
models On the square lattice, by using the saddle point equation we find 

Ngr = -2^, (36) 



which can be made more explicit by writing 



Ngr = 4vr^^ = Stt 

PsE{ps) 



l + -[\n- + 2]+0iz 



(37) 



where E is an elliptic function. 

All the above results can be expressed as functions of /? by solving the saddle point 
equation. Concerning asymptotic scaling, and therefore solving the saddle point equation 
at large P, one finds 

Mg ~ 4^2 exp (-y) • (38) 

The analytic structure of the various observables has been investigated in Ref. 0. The 
complex /3-singularities are square-root branch points, indeed quantities like x and behave 
as 



A{(3) + B{(3)^(3-(3, (39) 

around a singular point jSs, where A{(3) and B{(3) are regular in the neighborhood of f3s- 
The singularities closest to the origin are located ai (5 = 0.32162 (±1 ±i). Such singularities 
determine the convergence radius of the strong-coupling expansion, which is therefore (3r = 
0.45484, corresponding to a correlation length = 3.17160. 



C. The triangular lattice 

On the triangular lattice, using the results of App. ^ the saddle point equation can be 
written as 
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1 ^ r dki /-wvs 1 

t ~ i-TT 27r J-27T/V3 2-11 A{k) + z 



where 



1 COS ki + 2 cos — cos 

3 V 2 2 



(40) 



(41) 



and the momentum integration is performed over the Brillouin zone corresponding to a 
triangular lattice. By rather straightforward calculations (making also use of some of the 
formulas of Ref. p) the saddle point equation can be written as 



6 



-0 



1/4 



where 

Using the results of App. ^ one can find 



1 Z If 2\l/2' 



-1/2 



5 3z 3/, 2\i/2 

-H 1 + - 

2 8 2 V 6 



t 



X 



VgZ SPz 



-1/2 



(42) 



(43) 



(44) 



2 - 



1 

Z 



(45) 



u 



ml 



4(^+16 



(46) 



An estimator of the mass-gap fit can be extracted from the long distance behavior of the 
wall- wall correlation function defined in Eq. ([A7|) , indeed for a; 3> 1 



G^\x) oc e~ 



In the large- iV limit one finds 



M2 = ^ (^COsh^/it - 1^ 



z = Ml. 



A test of scaling is provided by the ratio 
l^t 2 



M, 



G 



"iz 



Arccosh 



3 

1 + :^^ 



^ 32^^ 10240' 



O^z") 



(47) 



(48) 



(49) 



where scaling violations are of the same order as those found on the square lattice for the 
corresponding quantity, cfr. Eq. (P7|). 



9 



The internal energy is given by the following expression 



1 ^ 
1 h -. 

6p 4' 



(50) 



leading again to the result (|32D for the condensate of the trace of the energy-momentum 
tensor, in agreement with universality. 

We calculated Qr on the triangular lattice, finding the following expression (in the deriva- 
tion we made use of the saddle point equation (HOl)) 



Ngr 



2 d\nz 



VS 0(3 ' 

which can be written in a more explicit form using Eq. 



(51) 



Ngr 



4:71 {1+-Z 



1/4 ^ 



E{pt) dpt 
1 - pi dz 



1 

24 



1 + -z 
6 



PtK{pt 



Stt 



1 + 2 (InA + 11) +0(^2, 
° ^ 48 6 ' ^ 



(52) 



where the continuum value of Ngr, obtained for 2; ^ 0, is in agreement with the results (|35[) 
and (13). 

In the weak coupling region t — the saddle point equation leads to the asymptotic 
scaling formula 



M, 



G 



4\/3exp ( 



2ti 



(53) 



The equations (|38| ) and (|53| ) are in agreement with the large- iV limit of the ratio of the 
A-parameters of the square and triangular lattice formulations calculated in App. |^, cfr. 
Eq. ( [A13| ), using perturbation theory. 

We have investigated the analytic structure in the complex /3-plane. Details of such 
study are presented in App. 0. As on the square lattice, the singularities are square-root 
branch points. Those closest to the origin are placed at /5 = 0.206711 ± i 0.181627, leading 
to a convergence radius for the strong-coupling expansion f3r = 0.275169, which corresponds 
to a correlation length C,g = 2.98925. 



D. The honeycomb lattice 

The analysis of models defined on the honeycomb lattice presents a few subtleties caused 
by fact that, unlike square and triangular lattices, not all sites are related by a translation, 
not allowing a straightforward definition of a Fourier transform. Nevertheless, observing 
that sites at even distance in the number of lattice links form triangular lattices, one can 
define a Fourier-like transformation that partially diagonalizes the Gaussian propagator (up 
to 2 X 2 matrices) [Q. In this section we present the relevant results, some details of their 
derivation are reported in App. 0. 
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Using the expression of i? ^ of Eq. (p^) we write the saddle point equation in the 
following form 



^ _ (3 _ fh dki f^/^ dk2 1 + 



where 



AW = g 



27r J-n/V^ 2ti A{k) + z (l + |z) 



2 - cos -^f^2 I COS -fci + COS -^f^2 



and integrating over the momentum we arrive at 

1/2 



1 _ P _ 1 

t 2^ 



1 + ^) Ph(^)ir(phW), 



where 



PhW = (l + ^j + (1 



(54) 



(55) 



(56) 



(57) 



From Eq. (|B4| ) we also derive 



X 



VsZ Zj3z 



(58) 



G 



(59) 



Z 

16 



(60) 



The two orthogonal wall-wall correlation functions G'^\x) and G^\x) defined in 
Eqs. (p6|) and (p^) allow one to define two estimators of the mass-gap from their long 
distance behavior 



Gf)(x) oce" 
G[7^(x) oce- 



(61) 



where x is the distance between the two walls in unit of the lattice spacing. In the continuum 
limit pv = fJ^h and they both reproduce the physical mass propagating in the fundamental 
channel. As on the square and triangular lattices, it is convenient to define the quantities 



Ml 



^ fcosh^ - 1 
9 V 2 

8 / , V3ph 
- cosh^-1 



(62) 
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which, in the continuum hmit, are also estimators of the mass gap. In the large- limit one 
finds 



Ml 



z. 



Notice that in the continuum large- limit the result 

M 



Mr 



(63) 



(64) 



where M is any mass-gap estimator, is found for all lattice formulations considered. 

On the honeycomb lattice the maximal violation of full rotational symmetry occurs for 
directions differing by a tt/G angle, and therefore, taking into account its discrete rotational 
symmetry, also by a 7r/2 angle. So a good test of rotation invariance of G[x) at large distance 
is provided by the ratio ^1^/11^: 



Arccosh 


1 + 1^(1 + 1 




-\/3Arccosh 


1 + 1^; 





640 



z' + Oiz- 



(65) 



As expected from the better rotational symmetry of the honeycomb lattice, rotation invari- 
ance is set earlier than for the square lattice, indeed the 0[z) scaling violation is absent. 
A test of scaling is provided by the ratio 



M, 



G 



32 



Arccosh 



1 9 

1 z H 

32 10240 



z 



+ O (^3) 



(66) 



where scaling violations are of the same order of those found on the square lattice for the 
corresponding quantity, cfr. Eq. (^). 
The internal energy is given by 



E 



1 z 
^-3^ + 4 



1 Ml 
1 \ 



(67) 



where the term proportional to Mq verifies again universality. 

In the weak coupling region t — the saddle point equation leads to the asymptotic 
scaling formula 



M, 



G 



4exp ( — — 



(68) 



The equations (pq ) and (pq ) are in agreement with the large-A^ limit of the ratio of the 
A-parameters of the square and triangular lattice formulations calculated in App. P, cfr. 
Eq. ( [B12| ), using perturbation theory. 

In Fig. 1^ we compare asymptotic scaling from the various lattices considered, plotting 
the ratio between Mq and the corresponding asymptotic formula (cfr. Eqs. (|38|), (^) and 
(|68|)). Notice that in the large- A^ limit corrections to asymptotic scaling are 0(Mg), in that 
corrections 0(1/ In M^) are suppressed by a factor 1/A^. 
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We have investigated the analytic structure in the complex temperature-plane of the 
N = oo model on the honeycomb lattice (details are reported in App. 0). As on the square 
and triangular lattices, singularities are square-root branch points, and those closest to the 
origin are placed on the imaginary axis at /3 = ±iO. 362095. The convergence radius for the 
strong-coupling expansion is associated to a quite small correlation length: C,g = 1-00002. 

III. CONTINUUM RESULTS FROM STRONG COUPLING 

A. Analysis of the series 

In this section we analyze the strong-coupling series of some of the physical quantities 
which can be extracted from the two-point fundamental Green's function. We especially 
consider dimensionless renormalization-group invariant ratios, whose value in the scaling 
region, i.e., their asymptotic value for f3 oo, concerns the continuum physics. Some strong- 
coupling series for selected values of are reported in the Apps. 0, ^ and |F| respectively 
for the square, triangular and honeycomb lattices. The series in the energy are obtained by 
inverting the strong-coupling series of the energy E = (3 + 0{(3^) and substituting into the 
original series in (3. 

Our analysis of the series of dimensionless renormalization group invariant ratios of phys- 
ical quantities, such as those defined in the previous section, is based on Pade approximant 
(PA) techniques. For a review on the resummation techniques cfr. Ref. 

PA's are expected to converge well to meromorphic analytic functions. More flexibility is 
achieved by applying the PA analysis to the logarithmic derivative (Dlog-PA analysis), and 
therefore enlarging the class of functions which can be reproduced to those having branch- 
point singularities. The accuracy and the convergence of the PA's depend on how well 
the function considered, or its logarithmic derivative, can be reproduced by a meromorphic 
analytic function, and may change when considering different representations of the same 
quantity. By comparing the results from different series representations of the same quantity 
one may check for possible systematic errors in the resummation procedure employed. 

In our analysis we constructed [//m] PA's and Dlog-PA's of both the series in (3 and in 
the energy. /,m are the orders of the polynomials respectively at the numerator and at the 
denominator of the ratio forming the [//m] PA of the series at hand, or of its logarithmic 
derivative (Dlog-PA). While [//m] PA's provide directly the quantity at hand, in a Dlog- 
PA analysis one gets a [//m] approximant by reconstructing the original quantity from the 
\l/m\ PA of its logarithmic derivative, i.e., a [Z/m] Dlog-PA of the series A{x) = Y^^o^i^^ 
is obtained by 

rx 

A/m{x) = aoexp dx'D\ogi/^A{x'). (69) 

where Dlogi/j^A{x) indicates the [l/m] PA of the logarithmic derivative of A{x). 

We recall that a [l/m] PA uses n = I + m terms of the series, while a [l/m] Dlog-PA 
requires n = I + m + 1 terms. Continuum estimates are then obtained by evaluating the 
approximants of the energy series a.t E = 1, and those of the /3 series at a value of j3 
corresponding to a reasonably large correlation length. 
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As final estimates we take the average of tfie results from the quasi-diagonal (i.e., with 
/ ~ m) PA's using all available terms of the series. The errors we will display are just 
indicative, and give an idea of the spread of the results coming from different PA's. They 
are the square root of the variance around the estimate of the results, using also quasidi- 
agonal PA's constructed from shorter series. Such errors do not always provide a reliable 
estimate of the uncertainty, which may be underestimated especially when the structure of 
the function (or of its logarithmic derivative) is not well approximated by a meromorphic 
analytic function. In such cases a more reliable estimate of the real uncertainty should come 
from the comparison of results from the analysis of different series representing the same 
quantity, which in general are not expected to have the same analytic structure. 

In the following of this section we present the main results obtained from our strong- 
coupling analysis. Most of them will concern the = 3 case. 



B. The square lattice 

On the square lattice we have calculated the two-point Green's function up to 0(/?^^) 
from which we have extracted strong-coupling series of the quantities E, x, ^hi '^i 



already introduced in Sec. [11 B|, and of the ratios r = M^/M^, s = M^/Mq. Some of the 



above series for selected values of A^ are reported in App. |^. Our strong-coupling series 
represent a considerable extension of the 14th order calculations of Ref . [l^ , performed by 
means of a linked cluster expansion, which have been rielaborated and analyzed in Ref. |]Tl] . 



We also mention recent works where the linked cluster expansion technique has been further 
developed and calculations of series up to 18th order and 19th order |T^] for d=2,3,4 
have been announced. 

In order to investigate the analytic structure in the complex /5-plane we have performed 
a study of the singularities of the Dlog-PA's of the strong-coupling series of x and ^q. As 
expected by asymptotic freedom, no indication of the presence of a critical point at a finite 
real value of /3 emerges from the strong-coupling analysis of A^ > 3 models, confirming 
earlier strong-coupling studies |Tl|. The singularities closest to the origin, emerging from 
the Dlog-PA analysis of x cind ^g, are located at a pair of complex conjugate points, rather 
close to the real axis in the A^ = 3 case (where jS ~ 0.59±z0.16) and moving, when increasing 
A^, toward the N = oo limiting points (3 = 0.32162 (1 ± i). In Table |I| such singularities 
are reported for some values of A^. The singularity closest to the origin determines the 
convergence radius of the corresponding strong-coupling series. For example for A^ = 3 the 
strong-coupling convergence radius turns out to be Pr — 0.61, which corresponds to a quite 
large correlation length ^ ~ 65. We recall that the partition function on the square lattice 
has the symmetry f3 — > — /3, which must be also realized in the locations of its complex 
singularities. 

By rotation invariance the ratio r = M^/M^ should go to one in the continuum limit. 
Therefore the analysis of such ratio should be considered as a test of the procedure employed 
to estimate continuum physical quantities. In the large- A^ limit r = 1 at all values of (3. This 
is not anymore true at finite A^, where the strong-coupling series of and M| differ from 
each other, as shown in App. |^. From G{x) up to 0(/5^^) we could calculate the ratio r up 
to 0(/5^^). The results of our analysis of the series of r for A^ = 3 are summarized in Table 
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[T|. There we report the values of the PA's and Dlog-PA's of the i?-series at E = 1, and of 
those of the /^-series at f3 = 0.55, which corresponds to a reasonably large correlation length 
^ ~ 25. We considered PA's and Dlog-PA's with / + m > 11 and m > I > 5. The most 
precise determinations of r*, the value of r at the continuum limit, come from Dlog-PA's, 
whose final estimates are r* = 1.0000(12) from the i?-approximants, and r* = 1.0002(6) 
from the /3-approximants (at (3 = 0.55). The precision of these results is remarkable. 

For all > 3 the violation of rotation invariance in the large distance behavior of G{x), 
monitored by the ratio /is/z^d, turns out quantitatively very close to that at = cxd when 
considered as function of (in a plot the A^ = 3 curve of /ig/ /id versus as obtained from 
the strong-coupling analysis would be hardly distinguishable from the exact N = oo one), 
/is/ /id is one within about one per mille already at ^ ~ 4. 

Calculating a few more components of G{x) at larger orders (i.e., those involved by the 
wall- wall correlation function at distance 6 and 7 respectively up to 0(/3^^) and 0(/5^^)), we 
computed the ratio 



Ml 



5 = T72 (70) 



up to 0(/3^^), by applying the technique described in Refs. |T^ , [T5[| . We recall that at A^ = oo 
we found s = 1 independently of (5. No exact results are known about the continuum limit s* 
of the ratio s, except for its large- A^ limit: s* = 1. Both large- A^ and Monte Carlo estimates 



indicate a value very close to one. From a 1/A^ expansion [|T6|,|T7| 



^ 0.006450 ^ / 1 > 



Monte Carlo simulations at A^ = 3 p| gave s = 0.9988(16) at f3 = 1.7/3 = 0.5666... 
(^ ~ 35), and s = 0.9982(18) at p = 0.6 (^ ^ 65), leading to the estimate s* = 0.9985(12). 

In Table |^ we report, for A^ = 3, the values of PA's and Dlog-PA's of the energy and 
/? series of s respectively at ii^ = 1 and at jS = 0.55. We considered PA's and Dlog-PA's 
with I + m > 13 and m > I > 5. Combining PA and Dlog-PA results, our final estimates 
are s* = 0.998(3) from the ii^-approximants, and s* = 0.998(1) from the j3 approximants 
evaluated at (3 = 0.55, in full agreement with the estimates from the 1/A^ expansion and 
Monte Carlo simulations. With increasing A^, the central estimate of s* tends to be closer 
to one. 

The scaling- violation pattern of the quantity fig/MG for A^ = 3 is similar to the pattern 
for N = oo (cfr. Eq. (P7D), i.e., it is stable within a few per mille for ^ ^ 5. 

Another dimensionless renormalization-group invariant quantity we have considered is 
u = ml/ (x"^4)5 whose large- A^ limit has been calculated in the previous section, cfr. Eq. (pll). 
At finite A^ its continuum limit u* is not known. From the expression of the self-energy 



calculated up to 0(1/A^) [p!q, p!^ , p!^] , one can obtain 

(72) 



0.006198 ^ / 1 



. 1 
u = — 

4 A^ ' " 

It is interesting to notice that the 0(1/A^) correction in Eqs. (fflD and (0) is very small. 
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At = 3 the analysis of the 0(/3^^) strong-couphng series of u detected a simple pole 
close to the origin at Pq = —0.085545 for the /3-series, and at Eq = —0.086418 for the energy 
series, corresponding to Mq = —16.000, which, within the precision of our strong-coupling 
estimate, is also the location of the pole in the corresponding N = oo expression (0). Being 
a simple pole, this singularity can be perfectly reproduced by a standard PA analysis, and 
indeed we found PA's to be slightly more stable than Dlog-PA's in the analysis of u. The 
results concerning A^ = 3, reported in Table |V| (for PA's with / + m > 16 and m > I > 8), 
lead to the estimates u* = 0.2498(6) from the energy analysis, and u* = 0.2499(5) form the 
f3 analysis (at j3 = 0.55). The agreement with the large- A^ formula ( [7^ ) is satisfactory. In 
Fig. 0the curve u{E) as obtained from the [10/10] PA and the exact curve u{E) at N = oo 
(cfr. Eq. |2l]) are plotted, showing almost no differences. 

In Table [V^ we give a summary of the determinations of r*, s*, and u* from PA's and 
Dlog-PA's of the energy and /^-series. 

We mention that we also tried to analyze series in the variable 

lN/2iN(3) ^^^^ 



iN/2-i{N(3y 

which is the character coefficient of the fundamental representation. As for the iJ-series, 
the continuum limit should be reached at a finite value 2 — > 1, and estimates of r*,s* and 
u* may be obtained evaluating the approximants of the corresponding z-series at z = 1. 
We obtained results much less precise than those from the analysis of the i?-series. Maybe 
because of the thermodynamical meaning of the internal energy, resummations by PA's and 
Dlog-PA's of the i?-series turn out much more effective, providing rather precise results even 
at the continuum limit E = 1. 

The strong-coupling approach turns out to be less effective to the purpose of checking 
asymptotic scaling. In Table [VI1| , we compare, for A^ = 3,4,8, as obtained from the 
plain 21st order series of and from its Dlog-PA's with some Monte Carlo results available 



in the literature. Resummation by integral approximants |20| provides results substantially 
equivalent to those of Dlog-PA's. For A^ = 3 Dlog-PA's follow Monte Carlo data reasonably 
well up to about the convergence radius [3r — 0.6 of the strong-coupling expansion, but they 
fail beyond I3r- On the other hand it is well known that for A^ = 3 the asymptotic scaling 



regime is set at larger /5- values |^. More sophisticated analysis can be found in Refs. |p!l| , |22| , 
but they do not seem to lead to a conclusive result about the asymptotic freedom prediction 
in the 0(3) a model. At larger A^, the convergence radius decreases, but on the other hand 
the asymptotic scaling regime should be reached earlier. At A^ = 4 and N = 8 the 21st order 
plain series of provides already quite good estimates of within the convergence radius 
when compared with Monte Carlo results. Again Pade-type resummation fails for j3 > j3r- 
We mention that at A^ = 4 the convergence radius (3r — 0.60 corresponds to ^g — 25, and 
aX N = 8 (3r — 0.55 corresponds to ^g — 8. 

In order to check asymptotic scaling we consider the ratio As/A2;, where Ag is the effective 
A-parameter which can be extracted by 
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where M is an estimator of the mass-gap, Rg is the mass- A parameter ratio in the square 
lattice nearest-neighbor formulation [ p3[] 



A' 



MS 



-2 



X v32exp 



As/ \ej r(l + ^ 
and A21 is the corresponding two-loop formula 



, 27tN \— f 2ttN ^ 
^21 = ( — — 7:P ) exp ( — -P 



2{N -2) 



(75) 



(76) 



The ratio As/A2z should go to one in the continuum limit, according to asymptotic scaling. 
The available series of Mq are longer than any series of the mass-gap estimators; therefore, 
neglecting the very small difference between Mq and M (we have seen that for > 3 
[Mq — M)/M ^ 10~^ in the continuum limit), for which formula (^) holds, we use Mq 
as estimator of M. In Fig. ^ we plot Kg/ ^21 for various values of N, N = 3,4,8, and for 
comparison the exact curve for N = 00. As already noted in Ref. by a Monte Carlo 
study, for A^ = 3, 4 at ^ ~ 10 the asymptotic scaling regime is still far (about 50% off at 
A^ = 3 and 15% at A^ = 4), while for A^ = 8 it is verified for ^ 4 within a few per cent, and 
notice that the convergence radius (3r — 0.55 corresponds to ^ ~ 8. Anyway with increasing 
A^ curves of As/A2/ clearly approach the exact N = 00 limit. 



C. The triangular lattice 

On the triangular lattice we have calculated the two-point Green's function up to 0(/?^^), 
from which we have extracted strong-coupling series of the quantities E, x, ^q, u, M^, 
already introduced in Sec. |11 (J| , and of the ratios s = M^/Mq. Some of the above series for 



A^ = 3 are reported in App. |^. 

Like 0{N) a models on the square lattice, no indication of the presence of a critical 
point at a finite real value of (3 emerges from the strong-coupling analysis for A^ > 3. By 
a Dlog-PA analysis of the 0(/5^^) strong-coupling series of x ^"^^ 'Cg at A^ = 3, we found 
that the singularities closest to the origin is /5 ~ 0.358 ± iO. 085, giving rise to a convergence 
radius f3r — 0.37 which should correspond to a rather large correlation length: C,g — 70. 
Increasing A^ such singularities move toward their N = 00 limit f3 = 0.206711 ± i 0.181627. 
Some details of this analysis are given in Table 

In our analysis of dimensionless quantities we considered, as on the square lattice, both 
the series in the energy and in f3. The estimates concerning the continuum limit are obtained 
by evaluating the approximants of the energy series ai E = 1, and those of the /9-series at a 
P associated to a reasonably large correlation length. For A^ = 3 we chose (3 = 0.33, whose 
corresponding correlation length should be ^ ~ 22, according to a strong-coupling estimate. 

Calculating a few more components of G{x) at larger orders (i.e., those involved by the 
wall- wall correlation function at distance \/3/2 x 5 up to 0(/?^^)), we computed the ratio 
s = /Mq up to 0{(3^^) [|1^,0. For A^ = 3 the analysis of the strong-coupling series of s 



some details are given in Table [VIII|) leads to the estimate s* = 0.998(3) from the energy 



approach, and s* = 0.998(1) evaluating the approximants at /3 = 0.33 (we considered PA's 
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and Dlog-PA's with / + m > 8 and m > I > 4). Such results are in perfect agreement with 
those found for the square lattice. 

PA's and Dlog-PA's (with / + m > 11 and m > / > 5) of the strong-coupling series of u 
expressed in terms of the energy, evaluated at E = 1, lead to the estimate u* = 0.249(1) at 

= 3. The analysis of the series in (3 gives u* = 0.2502(4). Again universality is satisfied. 

A summary of the results on the triangular lattice can be found in Table 0. 

As on the square lattice we checked asymptotic scaling by looking at the ratio Kt/ ^2h 
where At is the effective A-parameter on the triangular lattice, defined in analogy with 
Eq. ([7^. Beside the formulas concerning asymptotic scaling given for the square lattice 
case, cfr. Eqs. ([r^-|76|), we need here the A-parameter ratio At /As calculated in App. ^ cfr. 
Eq. (|A13|) . We again used Mq as approximate estimator of the mass-gap M. Fig. | shows 
curves of Kt/^2i for various values of N, N = 3,4,8, and for comparison the exact curve 
for N = oo. Such results are similar to those found on the square lattice: for A^ = 3,4 
asymptotic scaling regime is still far at C,g — 10, but it is verified within a few per cent 
at A^ = 8, where the correlation length corresponding to the strong-coupling convergence 
radius is ^ ^8. 



D. The honeycomb lattice 

On the honeycomb lattice we have calculated the two-point Green's function up to 
0(/5^°), from which we extracted strong-coupling series of the quantities E, x, ^q, u, M^, 
Ml, already introduced in Sec. |IID| , and of the ratios r = M^/M^, s = M^/Mq. Some of 
the above series for A^ = 3 are reported in App. ^ 

At A^ = 3 a Dlog-PA analysis of the 0(/3^°) strong-coupling series of x ^"^^ detected 
two couples of complex conjugate singularities, one on the imaginary axis at /5 ~ ±i0.460, 
quite close to the origin, and the other at /3 ~ 0.93±z0.29. The singularity on the imaginary 
axis leads to a rather small convergence radius in terms of correlation length, indeed at /? ~ 
0.46 we estimate i ~ 2.6. At A^ = 4 we found (3 ~ ±i0.444, and (3 ~ 0.88 ± ^0.41. At larger 
A^ the singularities closest to the origin converge toward the N = oo value /3 = ±i 0.362095. 
Notice that, as on the square lattice, the partition function on the honeycomb lattice enjoys 
the symmetry j3 —>■ —j3. 

Again we analyzed both the series in the energy and in (3. The estimates concerning the 
continuum limit are obtained by evaluating the approximants of the energy series aX E = 1, 
and those of the /3-series at /5 = 0.85 for the A^ = 3 case, which should correspond to ,^ ~ 22. 

By rotation invariance the ratio r = M^/M^ should go to one in the continuum limit. 
From G{x) up to 0{(3^^) we extracted the ratio r up to 0(/3^°). Again PA's and Dlog-PA's 
of the energy series evaluated at E = 1 and of the /3-series evaluated at /? = 0.85 (some 
details are given in Table ^) give the correct result in the continuum limit: respectively 
r* = 1.00(1) and r* = 1.001(1) at A^ = 3 (we considered PA's and Dlog-PA's with l + m > 16 
and m > I > 7). 

Calculating a few more components of G{x) at larger orders (i.e., those involved by 
G^\x), defined in Eq. (P7|), at distances x = \/3/2 x 9 and x = \/3/2 x 10 respectively at 
0(/334) and 0{P^^)), we computed the ratio s = M^/M^ up to 0{p^^) For A^ = 3 



the analysis of the strong-coupling series of s gives s* = 0.999(3) from the i?-approximants 
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and s* = 0.9987(5) from the /3-approximants evaluated at /3 = 0.85 (some details are given 
in Table XI), in agreement with the result found on the other lattices. We considered PA's 
and Dlog-Pa's with I + m > 22, m > I > 10. 

The analysis of the energy series of u confirms universality: PA's and Dlog-PA's (with 
I < m, I + m > 26, / > 12) of the energy series evaluated a.t E = 1 give u* = 0.249(3), and 
those of the /5-series at /3 = 0.85 u* = 0.2491(3). As for the square lattice, the curve u{E) 
obtained from the PA's at = 3 and the exact curve u{E) at N = oo, cfr. Eq. (0), would 
be hardly distinguishable if plotted together. 

As noted above, the convergence radius is small in terms of correlation length for 
all values of A^: it goes from ,^ ~ 1.0 at A^ = oo to ^ 2.6 at A^ = 3. Nevertheless in 
this case Dlog-PA resummations seem to give reasonable estimates of even beyond (3r 
(apparently up to about the next singularity closest to the origin). In Fig. ^ we show curves 
of Ah/A2i, where A^ is the effective A-parameter on the honeycomb lattice, for various values 
of N, N = 3,4,8, and for comparison the exact curve for A^ = oo. The necessary ratio of 
A-parameters has been calculated in App. p, cfr. Eqs. (|B11| ) and (pi2| ). 



E. Conclusions 

We have shown that quite accurate continuum limit estimates of dimensionless 
renormalizat ion-group invariant quantities, such as s and u (cfr. Eqs. (^) and (pi])), can 
be obtained by analyzing their strong-coupling series and applying resummation techniques 
both in the inverse temperature variable (3 and in the energy variable E. In particular, in 
order to get continuum estimates from the analysis of the energy series, we evaluated the 
corresponding PA's and Dlog-PA's at ii^ = 1, i.e., at the continuum limit. This idea was 
already applied to the calculation of the continuum limit of the zero-momentum four-point 
coupling gr, obtaining accurate results 0. These results look very promising in view of a 
possible application of such strong-coupling analysis to four-dimensional gauge theories. 



The summary in Table of our A^ = 3 strong-coupling results for the continuum 
values r*, s* and u*, for all lattices we have considered, shows that universality is verified 
within a precision of few per mille, leading to the final estimates s* ~ 0.9985 and u* ~ 
0.2495 with an uncertainty of about one per mille. The comparison with the exact N = oo 
results, s* = 1 and u* = 1/4, shows that quantities like s* and u*, which describe the 
small momentum universal behavior of G{p) in the continuum limit, change very little and 
apparently monotonically from A^ = 3 to A^ = oo, suggesting that at A^ = 3 G{p) is 
essentially Gaussian at small momentum. 

Let us make this statement more precise. In the critical region one can expand the 
dimensionless renormalization-group invariant function 

L(//M^) = £M (77) 
G{p) 

around y = /Mq = 0, writing 

Liy) = l+y + l{y) 

oo 

liy) = Y^c,y\ (78) 

1=2 
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l{y) parameterizes the difference from a generalized Gaussian propagator. One can easily re- 
late the coefficients Cj of the expansion ([78|) to dimensionless renormalization-group invariant 
ratios involving the moments of G{x). 
It is worth observing that 

u* = — ^ -. (79) 

4(1 -C2) ^ ^ 

In the large- limit the function l{y) is depressed by a factor 1/N. Moreover the coefficients 
of its low-momentum expansion are very small. They can be derived from the expansion 



of the self-energy [0,|T^,|T^ . In the leading order in the expansion one finds 

0.00619816... 



C2 

C3 
C4 



N 

0.00023845... 

N 

0.00001344 



5 



N 

0.00000090... , , 

C5 ^ , (80) 

etc.. For sufficiently large N we then expect 

Q < C2 < 1 for i> 3. (81) 

As a consequence, since the zero of L{y) closest to the origin is yo = —s*, the value of s* is 
substantially fixed by the term proportional to (p^)^ in the inverse propagator, through the 
approximate relation 

s* - 1 ^ C2 ~ 4m* - 1. (82) 

Indeed in the large- limit one finds from Eqs. (^) and (|72D 

, ^ , -0.000252 ^ / 1 \ 

where the coefficient of the 1/A^ term is much smaller than those of s* and u*. 

From this large- A^ analysis one expects that even at A^ = 3 the function l{y) be small in a 
relatively large region around y = 0. This is confirmed by the strong-coupling estimate of u*, 
which, using Eq. ([75|), leads to C2 — —0.002. Furthermore, the comparison of the estimates 
of s* and u* shows that s* — Au* ~ within the precision of our analysis, consistently with 
Eq. (^Tj). It is interesting to note that similar results have been obtained for the models 
with N < 2, and in particular for the Ising Model, i.e., for A^ = 1, where the strong-coupling 
analysis turns out to be very precise 0. 

We can conclude that the two-point Green function for all A^ > 3 is almost Gaussian in a 
large region around = 0, i.e., |p^/M^| ^ 1, and the small corrections to Gaussian behavior 
are essentially determined by the (p^)^ term in the expansion of the inverse propagator. 

Differences from Gaussian behavior will become important at sufficiently large momenta, 
as predicted by simple weak coupling calculations supplemented by a renormalization group 
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resummation. Indeed the asymptotic behavior of G{x) for x <^ 1/M (where M is the mass 
gap) turns out to be 

bo and 71 are the first coefficients of the /3-function and of the anomalous dimension of the 
fundamental field s respectively. Let us remind that a free Gaussian Green's function be- 
haves like ln(l/a:;). Important differences are present in other Green's functions even at small 
momentum, as shown in the analysis of the four-point zero-momentum renormalized cou- 
pling, whose definition involves the zero-momentum four-point correlation function ( p4D . 
However monotonicity in seems to be a persistent feature. 

Our strong-coupling calculations allow also a check of asymptotic scaling for a relatively 
large range of correlation lengths. For all lattices considered the ratio between the effective 
A-parameter extracted from the mass-gap and its two-loop approximation, A/A21, when 
considered as function of C,g, shows similar patterns with changing A^. Confirming earlier 
Monte Carlo studies, large discrepancies from asymptotic scaling are observed for A^ = 3 in 
the range of correlation lengths we could reliably investigate, i.e., ^ ^ 50. At A^ = 8 and 
for all lattices considered, asymptotic scaling within a few per cent is verified for ^ ^ 4, and 
increasing A^ the ratio A/A21 approaches smoothly its A^ = 00 limit. 
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APPENDIX A: THE TRIANGULAR LATTICE 

The sites x of a finite periodic triangular lattice can be represented in Cartesian coordi- 
nates by 

xih.h) = hVi + km, 
/i = 1, ...Li, I2 = 1, ■■■L2, 

r/i = (l,0), V2=(^l,^y (Al) 

We set a = 1, where the lattice space a is the length of a link. The total number of sites, 
links, and triangles is respectively A^^ = L1L2, Ni = SNg, Nt^ = ^N^. Taking into account 
periodic boundary conditions, a finite lattice Fourier transform can be defined by 

X 

k 

where Vg = \/3/2 is the volume per site, and the set of momenta is 
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r, , 27r ^ 27r ^ 

k{mi,m2) = —niipi + —m2p2, 
Li L2 

mi = 1, ...Li, 7712 = 1; ■■■L 



Pi 



P2 



2, 



0, 



Notice that 



Li L2 



(A3) 



(A4) 



To begin with, let us discuss the Gaussian model on the triangular lattice, which is 
defined by the action 



(A5) 



links 



where xi, Xr indicate the sites at the ends of each link. Performing the Fourier transform 
one derives the propagator 



AiKl =4 1 COS /ci + 2 COS — COS 

^ ' 3 V 2 2 



(A6) 



From these formulae one can easily obtain the large- limit, since the model becomes 
Gaussian for N 00. 

In the massive Gaussian model one may define an exponentiated wall-wall correlation 
function by 



(A7) 



In the O(A^) a models, in order to evaluate the ratio between the A-parameter of the 
MS renormalization scheme and the triangular nearest-neighbor lattice regularization, we 
calculated the correlation function G{x) in perturbation theory. In the x-space we obtained 
(neglecting O(a^) terms) 



G{x) = 
F{a/x) 



1 + ^^^tF{a/x) + O(t') 



N 



1 

2^ 



ln^-7,-ln2-lln3). 



(A8) 



In the p-space 



G{k) = 
D{ak) 



N -1 t 



N fc2 
1 



1 + 



A^ 



D{ak) 



2^3, 



27T 



Inafc -21n2 - -ln3 
2 



(A9) 
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The above results required the calculation of the following integral 

' F{a/x) + 0{a/x), 

2D{aq) + 0{aq), 



-TT 2tx J-2n/V3 27r A(A;) 
^ dki /■WVs dk2 A(g) - A{k) - A{k + q) 



TT 2tX J-2tt/V^ 271 



A{k)A{k + q) 



(AlO) 



where the extremes of integration are chosen to cover the appropriate Brillouin zone, which 
can be determined from the finite lattice momenta (|A3|). 

Comparing the above results with the two-point Green's function renormalized in the 
MS scheme 



G- 



G- 



MS V"^^ 

'k_ 



2e-^E 



MS 



1 



N 



- 1 t 



1 + 0{t 



N k^ 

and following a standard procedure, one can determine the ratio of the A-parameters 

1 \ 



A, 



4^3 



7r 



exp 



XN-2) 



(All) 



(A12) 



where At is the A-parameter of the 0{N) a models on the triangular lattice. Furthermore 
comparing with the Green function calculated on the square lattice P3[ one can also derive 



As 
At 



exp 



7r 



N-2 



V3 



r 



(A13) 



where As is the A-parameter on the square lattice. 



APPENDIX B: THE HONEYCOMB LATTICE 

The sites a; of a finite periodic honeycomb lattice can be represented in Cartesian coor- 
dinates by 

X = x' + pfjp 

X' = hffi + l2f]2, 

h = 1, ...Li, I2 = 1, ...L2, p = 0, 1, 

m=(^,^), V2={0,V3), v, = il,0). (Bl) 

We set a = 1, where the lattice space a is the length of a link. The total number of sites, 
links and hexagons is respectively Ng = 2L1L2, Ni = 3L1L2, Ah = L1L2. The coordinate p 
can be interpreted as the parity of the corresponding lattice site: sites with the same parity 
are connected by an even number of links. 

The two sublattices identified by x+(/i, I2) = x{li, I2, 0) and x_(/i, I2) = I2, 1) forms 
a triangular lattice. Each link of the honeycomb lattice connects sites belonging to different 
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sublattices. Triangular lattices have a more symmetric structure, in that their sites are char- 
acterized by a group of translations. It is then convenient to rewrite a field = h,p) 
in terms of two new fields = and = defined respectively on the 

sublattices x+ and x_. Taking into account periodic boundary conditions, a finite lattice 
Fourier transform can be consistently defined M 



x± 



4>±{x±) 



VhLiL2 



-ik-x-i 



(B2) 



where = 3-\/3/2 is the volume of an hexagon, and the set of momenta is 

- 27r ^ 27r ^ 
k = — mipi + —m2p2 
Li L2 



mi = 1, ...Li, 1712 = 1, •••-^2 



Pi 



^,0 



P2 



1 

'3' 



(B3) 



Using the results reported in Ref. 0] one can find the following expression for the matrix 
R-^ (cfr. Eq. (|13D): 



t 



E 



1 + 

Aki 



-ikj 



H{ky 



(B4) 



where 
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V3, f 3, 

2 — cos -— ^2 COS -fci + COS 



H{k) = e-''''^ (1 + 2e*t'=i cos ■ 



(B5) 



In the large- limit R~^{x;y) represents the two-point Green's function. 

In Ref. [14|, guided by the analysis of the Gaussian model on the honeycomb lattice, two 
wall-wall correlation functions were defined: 



G^:'\lh) = J2G{hfj, + i2m), 
h 

with the sum running over sites of positive parity forming a vertical line; 

4^^(173/) = Yl G{{1 - 2/2)771 + l2f]2 + Pffp), 

h,P 



(B6) 



(B7) 



where the sum is performed over all sites having the same coordinate X2- G^^{x) and G^\x) 
allows the definition of two estimators of the mass gap, fi^ and /^h, whose ratio must go to 
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one in the continuum limit by rotation invariance. On the honeycomb lattice the maximal 
violation of full rotational symmetry occurs for directions differing by a vr/G angle, and 
therefore, taking into account its discrete rotational symmetry, also by a 7r/2 angle. So a 
good test of rotation invariance is provided by the ratio between masses extracted from the 
long distance behaviors of a couple of orthogonal wall-wall correlation functions constructed 
with G{x) such as fiv/fJ'h- 

In the O(A^) a models, in order to evaluate the ratio between the A-parameter of the 
MS renormalization scheme and the honeycomb nearest-neighbor lattice regularization. we 
calculated, in perturbation theory, the correlation function 



Gix+ - y+) = (, 
In the x-space we obtained (neglecting O(a^) terms) 



(B8) 
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In the p-space 



G{k) 



N-1 t 
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l + -iD{ak) + 



3^3, 



D{ak) = — {\nak - 2 In 2) . 
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The relevant formulas required by the above calculations are can be found in Ref. 0]. 

Comparing the above results with the two-point Green's function renormalized in the 
MS, cfr. Eq. ( |A11| ), one can obtain the ratio of A-parameters 



A, 



4 exp 



TC 



N-2 3V3, 
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where Ah is the A-parameter of the honeycomb lattice, and also 

2 ^ 



— = V 2 exp 
A« 



7r 
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2 \ 2 



N-2\2 3VS, 

We also give few orders of the perturbative expansion of the internal energy 

- 1 t N -It^ 
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APPENDIX C: COMPLEX-TEMPERATURE SINGULARITIES AT AT = oo ON 
THE TRIANGULAR AND HONEYCOMB LATTICES 



In this Appendix we will compute the complex-temperature singularities for the N = oo 
model on the various lattices we considered. We will follow closely the analysis of Ref. 
for the square lattice. The N = oo solution is written in all cases in terms of a variable w 
related to the inverse temperature (3 by the gap equation which has the generic form 

(3 = h{w)p{w)K{p{w)) = b{w) (CI) 

for suitable analytic functions h{w) and p{w) and variable w, which will be defined below. 
Here K{w) is the complete elliptic integral of the first kind. In general b{w) will be defined 
on the complex plane with suitable cuts. To the purpose of computing the singularities 
of the inverse function w{j3) we will be forced to consider b{w) on its Riemann surface TZ. 
Moreover our discussion will be only local and thus we will determine all singularities which 
appear in the Riemann surface of w{f3). Notice that not all of them will necessarily appear 
in the principal sheet of w{P). 

Let us now consider a point Wq & TZ and let j3o = b{wo). To study w{P) in a neighborhood 
of Po let us expand b{w) around Wq. If b'{wo) is different from zero, then w{P) is obviously 
analytic in a neighborhood of (3q and admits an expansion of the form 

^ = ^0 + - l^o) + OiiP - (3o)'). (C2) 

Instead if b'{wo) = /3o is a singular point. Indeed let k be the smallest integer such that 
6W(wo) ^ 0. Then in a neighborhood of Po we have 

w=Wo+ [u^A - PoY" + - Pof"). (C3) 

and therefore /?o is a fcth-root singular point of w{f3). Thus in order to determine the 
singularities of w{b) we must determine the zeroes of V{w) on the Riemann surface of the 
function b{w). 

Beside the expansion in (3 we will be interested in expanding our observables in terms 
of the energy E. We want thus to study the singularities of the various observables when 
expressed in terms oi E. In practice we must study the functions w{E) and (3{E). 

For all lattices we consider we have 

E = -. + e{w) = rfr + e{w) (C4) 

where a is a constant and e{w) an analytic function. We will verify in each specific case 
that the zeroes of b{w) do not give rise to any singularity. Then by an argument completely 
analogous to the one given for w{(3), the singularities of 'w{E) are determined by the zeroes 
of ^ over the Riemann surface of the function E(w), which, because of the fact that e{w) 
is a simple rational function of w, coincides with the Riemann surface of b{w). 

Finally let us discuss P{E). Of course locally we can rewrite it as (3{w{E)). We shall 
show that dj3/dw and dE/dw never vanish at the same point. Then it is simple to convince 
oneself that the singularities of P{E) coincide with those of w{E). 
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1. Analytic structure of the complete elliptic integral K{w) 



The analytic properties of the function K{w) are well known (cfr. e.g. Refs. 
First of all K{w) is really a function of w^; indeed it has a representation in terms of 



hypergeometric functions as |8-27 



K{w) = ^2F,(^\,\;l;w'y (C5) 

From this representation one may see that K{w) is analytic in the w^-plane cut along the 
real axis from one to infinity. We want now to discuss the extension of K{w) to its Riemann 
surface. First of all let us introduce the theta function 

oo 

e3{v\T) = e"("''"+2'"^\ (C6) 

m=—oo 

63{v\t) is an entire function of r and v for Imr > 0. Moreover it satisfies the following 
properties (cfr. Ref. Chapter 5) 



es{v\T + 2) = e3{v\T), (C7) 

93{v/t\ - 1/r) = v^e™'/" ^3(t^|r). (C8) 



In terms of 63{v\t) we define the modular function 

^(") = ^ ^|(0|r) 

which is also an entire function of r for Imr > 0. The function A(r) has several important 
properties (cfr. Ref. ^6], Chapter 7): 



(i) Consider the group r2 of transformations 

ar + b 
CT + d 



r' = ^ (CIO) 



with 

M 



The function A(r) is invariant under i.e., A(r) = A(r'). Let us notice that r2 is generated 
by the transformations 

r' = r + 2 (C12) 

t' = . (C13) 

2r + 1 ^ ^ 

(ii) Let D be the domain of the complex plane bounded by the lines Rer = ±1 and the 
circles |r ± 1/2| = 1/2, including the boundaries with Rer < 0. Then for every r' with 
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Imr' > there exists a unique t E D and a transformation (|C10| ) connecting r and r'. Thus 
D is the fundamental domain of the group 

(iii) If c is a complex number different from and 1, the equation A(r) = c has one and 
only one solution in D. Moreover A(x + ioo) = for real x. 

Using the function A(r) we obtain a complete parameterization of the Riemann surface 
of K{w). Indeed, because of the last property, we can perform the change of variables 

= A(r), r G -D; then it can be shown that 



= lel{0\r). (C14) 

The complete Riemann surface is then obtained by letting r vary in the whole upper complex 
plane. Notice that since ^3(0|r) never vanishes for Imr > 0, the analytic extension of K{w) 
is always non-zero. 

We want now to express the analytic extension of K{w) on every Riemann sheet in 
terms of functions defined on the principal sheet. This will allow us to go back using w"^ as 
fundamental variable. This is easily accomplished using the second property of A(r). Indeed 
if t' is a generic point with positive imaginary part, there exists a unique t & D such that 
( piO| ) holds for suitable integers a, 6, c, d satisfying ( |C11| ). Thus 




mr') =el{Q\'^-^\. (C15) 



To compute the r.h.s., let us notice that the transformations ( |(J1CI| ) are generated by ( |U12[ ) 
and ( |C13|) . Then, using ( |C7| ) and 



which can be derived from (|C^ ) and ([C^), we get easily, by induction, 

^^3^(0|rO = . (demr) + ^0^3^(01 - 1/r)) (C17) 

where s assumes the values ±1. The presence of this sign is due to the fact that the matrix 
in ( |C11|) which corresponds to the transformation (|C10|) is defined up to a sign. If we define 



and fix the signs in M so that M is a product of Ti and T2 and their inverses only, then 
s = 1. Finally, using A(— 1/r) = 1 — A(r) which easily follows from (p8[), we get 



TT 



-^2^0|r') = s (dK{w) + icK{^l-w'^)) (C19) 

with w"^ = A(r'). As we shall see the sign s plays no role in the subsequent discussion. 

We have thus reached the following result J^: the Riemann surface of K{w) is obtained 
by considering matrices M generating transformations belonging to as M is defined 
modulo a sign we can always assume d > 0. Then each sheet of the Riemann surface is 
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labeled by the pair {d, c) of one row. (We write {d, c) instead of (c, d) to use the same 
notation of Ref. 0.) The extension K(^d,c)(yU!) on this sheet is given by 

s{dK{w) + tcK{Vl-w^)) (C20) 

where s is a sign depending on M. 

Finally let us notice that K{y/z*) = K{^/z)* so that K^cl,c){Vz*) = K(^d~c){\^)* ■ This 
observation will allow us to consider only the case d,c> 0. 



2. Location of the singularities 



The singularities of w{P) on the square lattice have been studied in Ref. |0. Here we 
will restrict our attention to the case of the triangular and honeycomb lattice. Let us first 
simplify the expressions (^2|) and (|5^) introducing new variables 



w 



w 



6, 
4 



on the triangular lattice 
on the honeycomb lattice. 



Then 



/3- 

where for both lattices 



1= pKip) on the triangular lattice 

— a/Sw pK{p) on the honeycomb lattice 
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(3u;-l)3/2(«; + 1)1/2 • 
The gap equation has an important property: as 

p{wy 



(C21) 
(C22) 

(C23) 
(C24) 

(C25) 



pi-wY 



p{wY — 1 



(C26) 



and the elliptic integral satisfies the property K{iz/ \/l — z^) = \/l — z'^K{z) we can im- 
mediately derive that (3{—w) = ±P{w) where the upper (resp. lower) sign refers to the 
triangular (resp. honeycomb) lattice. 

Let us now discuss the Riemann surface of b{w). As p^ is a meromorphic function with 
two poles in the w-plane for w = 1/2 and w = —1, we can apply the discussion of the previous 
paragraph. We must then consider the prefactor in front, which contains a square- root with 
two branching points and which has thus a double-sheeted Riemann surface. The Riemann 
surface is then labeled by two integers {d, c) as we discussed in the previous paragraph, and 
a sign s which specifies the sheet of the Riemann surface of the prefactor. Thus we have 
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1 1 



's,d,c 



zpK(^d,c){p) on the tr. latt. 



27T V3w 

f3s,d,c = s—V3wpK(^d,c){p) on the honeycomb lattice 
To identify the singularities we must then find the zeroes in the complex plane of 



's,d,c 



sp 



dw 



Attw (3w)±i/2 



±^(d,c)(p) + E(^d,c){P) 



(C27) 
(C28) 

(C29) 



where the upper (resp. lower) sign refers to the triangular (resp. honeycomb) lattice and 



E. 



(d,c) 



dE{w)+ic K{Vl 



(C30) 



The zeroes of ( |C29| ) have been studied numerically as in Ref. 0]. In Table pClll| we report 



the solutions with positive real and imaginary part we have found for the lowest values 
of ((i, c). We have verified that in all cases d'^l3s^d,c/ dw"^ 7^ at the singularity: thus all 
points are square-root branch points. Notice that if /3 is a singularity also — /5 and ±/3* 
are singularities. Our results are somewhat different from those of Ref. |0] on the square 
lattice. Indeed in our case we have found more than one zero with Re /? > and Im /3 > on 
each sheet. The principal sheet is an exception, as it is free of singularity for the triangular 
lattice, while it contains a pair of purely imaginary singularities for the honeycomb lattice. 
We stress that our search for zeroes of (|C29| ) has been done numerically and thus we cannot 
exclude that some zeroes have been overlooked. However we are confident that at least in 
the region < 2 our list is exhaustive for the values of (d, c) we have examined. 

To check the value of the zero with lowest we can compare with a direct determination 
of the singularity from an analysis of the high-temperature series of x by Dlog-PA's or first- 
order inhomogeneous integral approximants (lA's). In Table pClV| we report the results of 
such an analysis. These numbers are in very good agreement with the exact results, although 
the spread of the Dlog-PA's largely underestimates the true error. This is probably due to 
the fact that Dlog-PA's are unable to reconstruct the exact singularity. Indeed for [3 —>■ Psing 
we have X = Xo + Xi{P ~ PsingY^"^ + as xo 7^ this behavior can be reproduced by lA's 
but not by Dlog-PA's. 

Let us also notice that for the triangular lattice only the singularities with positive real 
part appear in our analysis. 

With series with 30 terms or more it is also possible to get an estimate of a second 
singularity. For the triangular lattice the series with 30 terms has a second singularity at 
P = —0.698(6) ± iO. 776(5) which corresponds to the second point in the sheet with d = 1 
and c = 2. For the honeycomb lattice we get ±0.449(13) ± zO. 610(17) (30 terms) and 
±0.4504(7) ± iO. 5909(25) (60 terms). The lA's are less stable and the 30-term series do 
not yield any result. With 60 terms on the honeycomb lattice we get ±0.449(3) ± zO. 585(2) 
which is in perfect agreement with the exact result. 

Beside considering the series in /3 we have also considered series with the energy E as 
variable. In this case we must consider the zeroes of dE /dw. Explicitly for the three lattices 
we have (for the square lattice we take w = Ps): 
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dE 1 d(3 1 ^ , 

~r~ = ir^~, square latt. V^^l) 

dE 1 d/3 



dw 6/?2 dw 
dE _ 1 d(3 
dw 3/3^ dw 



+ 3w triangular latt. (C32) 
+ 1 honeycomb latt. (C33) 



It is evident from this formulae that dE /dw ^ where dp/dw = 0. Thus, as we said at the 
beginning of this Appendix, the analysis of dE /dw provides all singularities of w{E) and 
(3{E). We get for the nearest singularities: 

E = ±0.330261131671 i square latt. 

E = -0.290013856190 ± 0.138180553789 i triangular latt. 
E = ±0.303078379027 ± 0.402035415796 i honeycomb latt. 

For the square and honeycomb lattice the singularity appears on the principal sheet of E{w) 
while for the triangular lattice it belongs the the sheet with {d, c) = (1, 2). 

From the position of the singularities we can now compute the convergence radius of the 
high-temperature series on the real axis. In terms of the correlation length, we find that the 
series converge up to ^conv where: 

1. square lattice: /3-series: ^conv = 3.17160; i?-series: C,conv = 1.38403. 

2. triangular lattice: /3-series: ^conv = 2.98925; i?-series: ^conv = 1-66706. 

3. honeycomb lattice: /3-series: ^conv = 1.00002; _E-series: ^conv = 2.43450. 

The series converge thus in a very small /5-disk: however PA's and lA's are quite successful in 
providing good estimates in a larger domain of the /5-plane: indeed for the square (triangular, 
honeycomb respectively) lattice, series with 21 (15, 30 resp.) terms give estimates which 
differ from the exact result less than 1% till ^ ~ 10 (5,15 resp.). 



3. Conformal transformations 

Once the singularities are known one can use a conformal transformation to get rid of 
the nearest ones (cfr. e.g. 0) and thus accelerate the convergence of the approximants. 

Let us first consider the triangular lattice which has two singularities located at /? = pe^*^ 
with p = 0.27516911105, 6 = 0.720896055. As in Ref. we consider a transformation of 
the form 

I3 = pw[l + P{w)Q{w) + /i Qiwf] (C34) 

where Q{w) = 1 — 2wcos6 ± is a polynomial which vanishes for w = e^*^. P{w) is 
determined requiring that dj3/dw = for w = e^*^. The simplest polynomial with this 
property is given by 
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P(w) = Kr- (cos 26- w cos 6) . (C35) 

2 sin 9 



We will also consider a second possibility given by 
cos 46* cot 29 



(3 = pw 




) (1 -2«;^cos20 + w*) + ^(1 -2w^cos2^ + w 



sin 2^ 



(C36) 



This transformation would also work if we had four singularities at (3 = ±pe^*^. It is easy 
to see that for 6* = 7r/4 it reduces to the transformation used in for the square lattice. 
For the honeycomb lattice the nearest singularities are f3 = ±ip with p = 0.3620955333. We 
can thus use ( |C34| ) with 9 = 7r/2, i.e., 



P = pw 



(C37) 



In all cases p is a. free parameter which can be used to optimize the transformation. 

In order to compare the series with and without conformal transformation we have com- 
pared the results for the magnetic susceptibility for series with 15 and 30 terms on the 
triangular lattice and honeycomb lattice respectively. 

For the triangular lattice we have considered the [6/7], [7/6], [6/8], [7/7] and [8/6] DLog- 
PA's. For the series without conformal transformation we have found that all PA's have 
singularities on the real axis or with a small imaginary part {lm(3 ^ 0.2) with 0.4 ^ Re/5 ^ 
0.6. Excluding only the PA with the nearest singularity ([6/7] for which f3sing — 0.414) we 
find that the estimates agree within 1% till (3 ~ 0.33 (^ = 5.31) where we get x = 57.69(56) 
(this value is the average of the estimates of the various PA's while the error is the maximum 
difference between two different PA's) which must be compared with the exact value x = 
56.9837. Let us now consider the series obtained from the conformal transformation ( |C34| ) . 
The results are now much more stable: at /? = 0.33 the series with p = gives x = 56.978(25) 
while for p = 0.5 we get x = 56.9804(24). The estimates of the series with p = agree 
within 1% till /3 = 0.46 (^ = 21.58) where we get x = 672(6) (exact value x = 674.86) , 
while for the case p = 0.5 the same is true till (3 = 0.61 (^ = 110.27). In this last case 
however the fluctuations of the approximants are not a good estimate of the error: indeed 
we get X = 12861(126) to be compared with the exact value x = 13290. The estimates agree 
within 1% with the exact resuh only till p = 0.56 (^ = 64.00). 

Let us now consider the second transformation. The Dlog-PA's for p = do not have any 
singularity on the real axis for /3 < 3 and their estimates are extremely stable. For the value 
we have considered before, f3 = 0.33, they give x = 56.9808(5) which is in excellent agreement 
with the exact value although the error bar is clearly underestimated. The five different PA's 
agree within 1% till /3 = 0.72 (^ = 365.05) where the estimate is x = 110200(1100) to be 
compared with the exact value x = 123387. Again the error bar is underestimated and the 
true error is 11%. Indeed the estimate from the PA's agrees with the exact value within 1% 
only till /3 = 0.53 (^ = 46.18). 

If we use the series with p ^ the stability of the PA's decreases and singularities on the 
real axis begin to appear near the origin. For p = 0.1 the PA's are still reasonably stable 
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and wc get for (3 = 0.33, x = 57.03(6). PA's agree within 1% till p = 0.41 = 12.55) where 
we get X = 257.9(2.4) to be compared with the estimate with /i — 0, x — 255.83(2) and the 
exact value x = 256.02. 

An analogous analysis can be done for the honeycomb lattice. In this case we have 
considered the [14/14], [13/16], [14/15], [15/14], [16/13] Dlog-PA's excluding in each case 
those having a singularity on the positive real axis with f3 < 3. We find that PA's of the 
standard series agree within 1% till f3 = 1.15 = 16.23) where x = 303.2(2.8) to be 
compared with the exact result x = 305.53. The conformally transformed PA's give instead 
X — 305.14(69) for = and x — 305.84(63) for n — 0.3. The conformally transformed 
PA's agree within 1% till /3 = 1.35 (^ = 33.50) where we get x = 1100(15) for = and 
1114(12) for fi = 0.3. This should be compared with the exact value x — 1108.13 and with 
the estimate from the standard series x = 1073(47). 

From this analysis it emerges that the use of the conformal transformation is extremely 
useful, especially for the triangular lattice where one gets results which are better by a factor 
of 10 - 100. Of course the interesting problem would be to generalize the method to finite 
values of N. There are two problems here: first of all the exact location of the singularities is 
not known. Moreover also the exact nature of the singularity must be determined from the 
series. The first problem is probably not a very serious one: indeed if we redo the analysis 
we have presented using the values of the zeroes obtained from a Dlog-PA's analysis of the 
series itself the results are essentially unchanged. The really serious problem (at least for low 
values of N) is the nature of the singularity: indeed the transformations we have considered 
apply only to square-root branch points. If the singularity is different new transformations 
must be used. 



APPENDIX D: STRONG-COUPLING SERIES ON THE SQUARE LATTICE 

A complete presentation of our strong-coupling series is beyond the scope of the present 
paper. A forthcoming publication will include all the relevant "raw" series. In order to 
enable the interested readers to perform their own analysis, we present here just the series 
for the internal energy E. the magnetic susceptibility the three mass scales Mq, M^, and 
M^, and the renormalization-group invariant quantity u, on the square lattice, for the most 
interesting values of N, i.e., N — 3,4,8. The following Appendices will be devoted to the 
triangular and honeycomb lattice. 
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X = 1 + 4/5 + 12/52 + i§0/53 + 84/3^ + ^/^^ + if^/56 + i^/3^ + ^/38 + ^/?9 

I 213868 ^10 I 8483180 ,oll , 77643788 ql2 , 138079108 ,ol3 , 3569641036 ^14 , 17867938876 ql5 

27 I' 567 ^ 2835 ^ 2835 ^ 12525 ^ 127575 

I 1066237156 ^j16 , 71108!)8:-i028 ql7 , 992250411932 ^18 , 45911479386812 ql9 , 11971698881708 q20 
4725 ^ 212625 ~'~ 1913625 63149625 ^ ~'~ 12629925 

^ 69463206903148 ^21 ^Q(^22)^ (D2b) 

A/f2 _ /I _L 11 I 34 ;o3 4^4 224^5 80 /06 , 4 ;o7 544 ;o8 58150 /o9 , 75032 /0 10 

Mq- p - 4 + yP + -g-p - gP - Y35P - + ~ '45 y ~ TTOlf^ + ^05~^ 

_ 7465004^11 I 20575936^12 , 204232808 ql3 _ 359904856^14 , 34575314246 ^15 _ 938720152 ql6 

18225 P ~^ 42525 ^ ~'~ 382725 127575 ~'~ 5740875 382725 

_ 3116875126738 /ol7 , 47791394648 ^18 _ 128064328183586 o 19 , (n(q20\ (r>9r\ 
189448875 1148175 3978426375 ~l~ '^Kh' )i \U^(.) 

— /9-I _4_|_ll/3_|_M/33_4a4_3M«5_W^6_836^7|404^8^ 215318 q9 , 4700 qlO 
s ^"raA-'fg/-' a/-' 135/^ 3/^ 81/''^ "'" 405'"^ 8505 81 ^ 

_ 443038 ,011 I 24410738^12 _ 24732128 ,ol3 , 139715614^14 , 3172203256 ql5 , r)(qlG\ (D9r\'\ 
6075 127575 382725 ~'~ 1148175 ^ ~'~ 5740875 "r '-^Im J) l.-U^UJ 

/\^2 _ ^3-1 _ 4 I 11 ,Q I 65 /Q3 _ 3341 q5 _ 6781 q7 _ 1026007 q9 , 76323223 qll , 6716460083 /o 13 
iVJj fj 'I -r g /J -^ -^g/J 540 /-"^ 540 /-"^ 136080 ~'~ 2041200 36741600 

(D2e) 



+ 0(/3^^), 

, = 4/3 - 48^2 + if^/3^ - 6560/3^ + - 268|576^6 ^ 523™^7 _ 1833750976 ^8 

I 192827317268 q9 _ 2252963987728 qlO , 61420977655724 oil _ 2152899064881952 ql2 
' 135 135 315 945 

I 3593447617956212 yol3 _ 1469484644425188368 ql4 , 7358240125587547636 ql5 

135 4725 2025 /""^ 

_ 1805422873668489252736 glG , 316414144175527075263484 ^17 _ 11090798958461896523514928 ^18 

42525 637875 1913625 

_l_ 475138205142448133881786892 yol9 _ 49962959062246660123963410208 q20 
' 7016625 63149625 
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_|_ 583759338332900106172953943052 p21 _|_ q^012\^ (D2f) 



3. AT = 8 



_ fl_|_6|o3_38/Cj5_ 19924 q7 i 908/o9 , 103032376 /oll , 284005868 ql^ _ 801343307824 gib 

^ h> ^ ^h> x^h> 2625 "'" 525'^ 1771875 1771875 3410859375 

74654425999556 ,^17 5113078463377748 ,^19 , 1298678765126787388 ,^21 , (^( q22\ (r\'Xo,\ 

30697734375 1995352734375 69837345703125 "r "^VM ^ i^L/Od,; 

X = 1 + 4/3 + 12/3^ + ^/3^ + ^(3' + + ^/^^ + + 

_l_ 6597628 g^ , 160511644 yolO , 433554844^11 , 4869604484^12 , 17075166692^13 , 180200884868 14 

2625 1^ 39375 1^ 70875 1^ 590625 1771875 20671875 

_l_ 407707485308/^15 _ 18904178053196 gl& _ 861976694991428 on _ 2351633289032404 glS 
136434375 2046515625 30697734375 51162890625 

421417265877650188 ,:jl9 , n( (T^'XW\ 
9976763671875 P ^ '^KP ^ l^UODJ 
/Lf2 _ /O-l _4il9/0|298/o3_28/o4_ 5254 nb _ 2128 ^6 _ 654656 q7 , 98032 ^8 _ 95888554 q9 
IVlQ fJ ^ -|- J. -|- 25M 875 375 39375'^ ^ 5625 1771875'^ 

I 5421752^10 _ 1465956704 /Oil _ 2887857152 /0 12 , 21729605024744 ,313 _ 3782033297144 /ol4 

28125 P 20671875 P 8859375 ^ ' 10232578125 P 930234375 ^ 

I 354224260413706 /ol5 , 72788687225944 gie _ 204193924410409426 gl7 , 40434518800599992 0I8 

153488671875 P 7308984375 P 4655823046875 P 767443359375 P 

I 7755263425304938826 ,019 , nf^'^O') (Tt'^r') 
^ 241744658203125 P ~^ '^U'' > ■> \UO<^) 
T\yr2 _ 0-1 _4_|_19/0|298 ^3_28^4_ 5744 ,o5 _ 182 /06 _ 853316 ,0 7 , 343772 ,08 _ 25597438 ,o9 
s y ' 5 P ' 75 P 25 P 875 P 75 P 39375'^ ~'~ 28125'^ 1771875*^ 

I 10063732 A5l0 I 7958412578 /Oil , 637810442 al2 , 24367739089468 /o 13 _ 22983973787546 a 14 

140625 P ~^ 103359375 P ~'~ 44296875 P 51162890625 P 23255859375 P 

I 828277050985928 /315 , r^r q16\ (-n'iA\ 
767443359375 P + '^IP J' {U6a) 
/Vf2 I f)-l _4_|_19/0_|_23/o3_ 101219 /o5 _ 1598337 ,o7 , 3042153439 ,o9 , 716944244341 all 
IVI^ -t-/J ^-t-g/J-t-g/J -j^QgQQ /J g^gQQ /J -t- 141750000 4961250000 P 

I 668912806018181 ,013 , r)( gl5\ (D'^(^\ 
' 2728687500000 P ' J' {UO^) 

— AR /1S/^2 I 2804 ,33 32736 qA , 1146292 ^5 22299088 , 27328827644 q7 106327225472 q8 

+ —13 i^jJ + — 2625 — P 875 P 

1241049893108 ,o9 _ 72427584678352 ,olO , 48910909763139412 ,oll _ 190295816674721312 ,ol2 

875 P 4375 P 253125 P 84375 P 

46643731943589538628 ,ol3 _ 2117208265068678482704 glA , 83213961927741336432836 ,^15 
1771875 P 6890625 P 23203125 P 

_ 142777012206408305834115008 ,ol6 , 14998421484199502791600778428 ^17 
3110859.375 ■ ^ 30697731375 

2917K8910(il(i5!)0173820050!)93776 ol8 , 132815:-i95(i9200995091K19594751(i884 ,ol9 

51162890625 P ^ 1995352734375 P 

861233015742221409021216372394912 ,020 , 211098292286114752398400491787924004 ,o21 , n22\ (T)ot\ 
1108529296875 P 23279115234375 P ^\P )■ K'-''-''') 



APPENDIX E: STRONG-COUPLING SERIES ON THE TRIANGULAR 

LATTICE 

Wc present here the series for the internal energy the magnetic susceptibihty the 
two mass scales Mq and M'^, and the renormalization-group invariant quantity m, on the 
triangular lattice, for A?" = 3, 4, 8. 

1. iV = 3 

« _L O «2 , 17 ,o3 , /I ^4 , 132 ,o5 352 ,06 53833 ,o7 246072 gS 38879794 ,o9 6397316 ^jlO 

_ 8815166536 /Oil _ 63896967072 /ol2 _ 81729013664 /0 13 , 23710274924992 /0 14 , 415243711166548533 ,ol5 
2786875 P 11790625 P 15640625 P ' 2360483125 P ' 5016026640625 P 



U 
+ 
+ 
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+ C'(/3^^) (Ela) 

^=1 + 6P + 30/^2 + 6Z2^3 ^ 2|2^4 ^ 3^^5 ^ 14Z|Z84^6 ^ 3821432^7 ^ 68^^^8 

I 131380789212 ^9 , 449739783516 ^10 , 12494938177056 ^11 , 90716274919896 ^12 

336875 f ~'~ 336875 ^ ~'~ 2786875 ^ '~ 6131125 ^ 
I 36820597127739144 /ol3 , 9067857431725940688 ;Ol4 , 2430225926454897621936 /ol5 , /n//3l6\ fpiK^ 
' 766390625 I- ' 59012078125 ~'~ 5016026640625 I- ' ) \s^J-0) 

M"^ = — 4.-1-^/9-1-32/32, 4552 /o3 , 2176 oA _ 40312 /O 5 _ 288944 n& _ 13545688 o7 

IVlQ ^jj ^ -|- ^g/J -|- 5 /J ~r 525 525 2625 3675 67375 1^ 

_ 372621008 oS _ 338341556368 _ 33053110624 alO _ 5766016722184^11 , 279009107673392 ,ol2 

1010625 459834375 ^ 18393375 2299171875 59012078125 

I 72472337922926768 ,^13 , /n/'/3l4^ 

+ 2149725703125 ^ + > \^^^) 

M"^ = -R~^ —4-1- -I- ^/^^ -I- M96/Q3 I 11944 o4 13836 /o5 _ 6627032^6 _ 337042602 q7 

IVl^ g/J ^-t-^g/J-t-g/J "r 525 M "T 2625 875 1^ 91875 1684375 

_ 1079483686 /o8 _ 408531412013 a9 _ 3031149364472 glO , a^^\ fFIHl 
2358125 459834375 ^ 2299171875 ^ ) l^Hilu; 

K = 6/3 - 108/^2 + 9f2/33 _ 33840/3^ + IMgp/^s _ 371^^6 _^ 164533949256 ^7 

_ 20399809562064 a8 , 2838986498807112 q9 _ 50284702987803792 a 10 , 405247036983420442416 oil 

6125 48125 48125 21896875 '"^ 

_ 50244722315696297988096 ol2 , 4449723414146572825129368 /ol3 _ 866958410946002558379625488 /ol4 

153278125 766390625 8430296875 

_l_ 1827333381636135987619199322672 olb , /ni'/^^^N /pi \ 

1003205328125 '^\i-> ) \^^^) 



2. N = 4 



X = 
+ 
+ 



7? « -L O «2 I 10 ,03 , 10 ;o4 S a5 238 o6 4312 /q7 12316 ,38 30212 ,^9 60434 olO 

1579292 ,311 I 24363404 ,ql2 , 869521126 ^; 13 , 3674593642 ,3 14 , 34840110848 /0 15 , n( R^^\ CPOc,\ 
94^^ + 14175 P + 42525 P + 42525 P + 127575 P + ^\P ) \^^^) 

1 + 6/3 + 30/32 ^ ^34^3 ^ 554^4 ^ 2162/3^ + ^/36 + i^/3^ + 1^(3^ 

15035566 ,o9 , 3260594 olO , 361103026 all , 50061022798 ol2 , 450734165906 olZ 

15 Z P ^ 105 ^ 4725 P 14175 

439170936682 ^14 ^ 11246(^^5102 ^15 ^ ^^^16) ^E2b) 

- 4 + f /3 + f /32 + ^/33 + f /3^ - ^/35 - i||4^6 _ 3||24^7 

_ ^„J104 qS _ 19703848 q9 _ 45897388 qlO , 357801536 oil , 2167022524 /0 12 , 83057677936 /ol3 
405 P 25515 42525 382725 P 127575 P 1148175 

+ 0(/3i4) (E2c) 

A/f2 _ 2/3-1 _ /I I 34 ,o I 20 /o2 I 245 /o3 I 107 /o4 _ 18581 /35 _ 40037 gG _ 324218 ,^7 
^^H — zP 9P ^ ■iP ^ 27 P ^ 27 P %10 P 40b P 1215 P 

2772953 ,38 189804947 n9 493593379 ,310 , r^/ nll\ (-pnA\ 
- -^m^P - 204120 P - 765450 P + ^^P > \^^^) 

u^6p- 108(3^ + 1910/3^ - 33828/3^ + 599078/3^ - 10609380/3^ + 28i8307806 ^7 

_ 49910910532 ,38 , 2651696611738 /39 _ 46960305901996 ,310 , 17464545810431714 ,311 

15 P ^ 45 P 45 P ^ 945 P 

_ 220920665969769296 ,312 , 82160433528738212702 ,313 _ 1455022824164792634328 ,314 

675 P 14175 P 14175 P 

^ 5153554644036578704006 ^15 ^ ^^^16) (E2e) 



3. AT = 8 

p _ /3 I 9/32 I 10 /33 I 10 /34 8/35 238/36 4312 ,3? 12316 ,38 30212 ,39 60434/310 

h - /3 + 2/3 + y/3 + ^/3 - 3/3 - —/3 - -j^/3 - -^/3 - -^/3 - -^/3 

1579292 ,311 I 24363404 /312 , 869521126 ,313 1 3674593642 ,314 , 34840110848 ,315 , n( f^^^\ 
94^P + 14175 P + 42525 P + 42525 P + 127575 P + ^^/^ ) 1^"^^) 

X = 1 + 6/3 + 30/3^ + 134/3=^ + 554/3^ + 2162/3^ + ^/3« + i^/3^ + lM/38 

I 15035566 ,39 , 3260594 ,310 , 361103026 ,311 , 50064022798/312 , 450734165906 ,313 
45 P ^ 3 P ^ 105 P ^ 4725 P 14175 P 
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^ 439170936682 ^14 ^ 112466^18825102 ^15 ^ ^^^16) (gg^) 

]\/r2 _ 2 n~l A j_ 34 n j_ 20 n2 j_ 248 q3 j_ 32 n4 9268 /Q 5 14024 ^6 323204 ^7 198104 ^8 
^^^G - sP - 4 + -g-p + y/J + + yP - lo5"P ~ ^35~P ~ ^2T5~P 405~^ 

_ 19703848 o9 _ 45897388 ^10 , 357801536 ^11 , 2167022524 ^12 , 83057677936 ^13 , /n('/3l4\ ('pQp^ 

25515 P 42525 ^ ~'~ 382725 ^ ~'~ 127575 ^ ~'~ 1148175 ^ ~T~ J \i^OL) 
= "^R^^ — A 4- _L 20/o2 I 245/03 , 107/^4 _ 18581 q5 _ 40037 f>6 _ 324218 a7 _ 2772953 a8 

ivi^ g/j 'I -|- g -r g -r 27 P 27 P 810 P 405 ^ 1215 ^ 4860 ^ 

189804947 ^9 493593379 ^10 , ^/^;11n ff?OA\ 

- 204120 P - 765450 + ) l^"^^) 

M = 6/5 - 108/^2 + 1910/3^ - 33828/5^ + 599078/5^ - 10609380/5^ + 28i8307806 ^7 

_ 49910910532 ^8 , 2651696611738 o9 _ 46960305901996 ,olO , 17464545810431714 ^11 

15 P 45 P 45 y ' 945 P 

_ 220920665969769296 ;Ol2 , 82160433528738212702 ql3 _ 1455022824164792634328 ql4 

675 P ~^ 14175 P 14175 ^ 

^ 5153554644036578704006 ^15 + 0(^16) (£3^) 



APPENDIX F: STRONG-COUPLING SERIES ON THE HONEYCOMB LATTICE 

Wc present here the series for the internal energy E, the magnetic susceptibihty x, the 
three mass scales Mq, M^, and M^, and the renormalization-group invariant quantity u, on 
the honeycomb lattice, for N = 3,4, 8. 



1. N = 3 



P _ fQ_3o3:88o5_ 1761 o7 , 14902 q9 _ 3439769596 yoll , 2045078768 ,ol3 _ 5158924941321 ^15 
^ y 5/"^ ~'" 35^"^ 175 P ' 385 ^ 21896875 ^ ' 3128125 ^ 1861234375 ^ 



+ 



28899507590512836 ql7 _ 172552140904013898254 ,ol9 , 415753493847105514488 ^21 
2425932884375 ^ 3335657716015625 ^ ^ 1835409748046875 ^ 

_ 871596813626516225704857972 ,q23 , 297323210393498586435793176 q25 
872691149952587890625 ^ ~^ 67130111534814453125 ^ 
93896946096105935045256448288 ;o27 , 6327983006511909061486897651308648928 /o29 , r)( 

4757152640718994140625 ^ ' 71634557740672624383232421875 ^ ' '^KH ) ■'-^J 



X=l + 3P + 6(3' + f^' + f /?^ + f /3^ + ^/?6 + ^/5^ + '-^P 

I 52831068 ,010 I 1370148342 oil 1 36519432 ^12 , 2120228676 ^13 , 81696113448 ,ol4 

336875 P ' 21896875 ^ ' 398125 ^ ' 3128125 ^ ~^ 109484375 ^ 
_ 4005762392397^15 _ 5554706676084 am , 522250407043193142 ,ol7 , 28557036836995675104 ,ol8 _ 
2605728125 ^ 3648019375 ^ ' 60648322109375 ^ ' 3335657716015625 ^ 

84308531983195740042 ^19 _ 568341711070335425448 ,o20 , 60997222972929327707388 ,o21 , 
2382612654296875 ^ 16678288580078125 ^ 383600637341796875 P 



+ 



+ 



2382612654296875 ^ 16678288580078125 ^ ~'~ 383600637341796875 ^ ' 

26213984562109070694215784 q22 _ 32254495435296445139757114 n23 _ 43715635286204172906375312 ^24 

171538289990517578125 ^ 45931128911873016875 ^ 67130111534814453125 ^ 

546580 109 in55()ll:i()7() 102 1184 ,'j25 , 2495:-i7.S9:-i:i:-i 192()i:-ii:-S7.S2223032 ^26 

i715:-!8289990517578125 ' 872691119952587890(125 

682220109367777 129732001 181021372 n27 _ 72810 132()2112487309281>')1925S95608 o28 

48718000193603218994140625 ^ 57712092537037659423828125 ^ 

562809908853468932757594137982353385288 a29 , 8519985565279875808402224146574969010656 a30 
8954319717584078047904052734375 ^ 152223435198929326814368896484375 ^ 



+ 0{P^') (Fib) 

_ 4 A5-1 _4_|_52^_316^3_|_ 5472 nb _ §4 n6 _ 19122724 n7 , 47808 08 , 1316732616 q9 
IVlQ ^fj ^ -t- -^^fj ^^^fj -t- fj 1010625 P ~^ 4375 ^ ~'~ 21896875 ^ 

_ 7823616 olO _ 472105313248 ^11 , 1973349088 ^12 , 46661184606992 ^13 _ 26076544910176 q14 

153125 P 2299171875 ^ ~'~ 8421875 ^ ~'~ 65143203125 ^ 26823671875 P 

_ 1783696991213198452 ,ol5 , 3724145632458848 ql& , 5630519133419621984 nl7 _ 14678226956585492416 q18 

667131543203125 ^ ^ 938828515625 ^ ~^ 526682797265625 ^ 877804662109375 ^ 

_ 23307722455628565130807096 ql9 , 131635204541138381824 ,o20 , 14931213920583922857376688 q21 

523614869971552734375 ^ 1796123385546875 ^ 79335586359326171875 ^ 

_ 12490315953001109355570208 n22 _ 268173518152038112779584900384 ,o23 

37943106519677734375 ^ 335986208231746337890625 ^ 

_l_ 1303640079018735894344620352 ,^24 , 38043748942433272601810917870624 ^25 

872691449952587890625 ^ 11242615429293050537109375 ^ 

_ 35163452856807151015094953488864 q2& _ 32741212467941707889297214938863984773664 q27 
5174187606768893603515625 ^ 2283351527983939902215533447265625 ^ 
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1538444172762671916055304195595321312 ,028 , r>f/?29N (V^r) 
49801555715150600933837890625 ' '^KH J \^ ■'-^) 

2/Q-2 _ 38 _|_ 1418/32 , 2804^4 _ 4877702 n6 , 83484008 aS _ 510322072076 olO 
9*^ 45 ~'~ 1575^ ~'~ 7875 3031875^ ~'~ 21896875^ 34487578125 ^ 

241521726513358 ^12 35557059082336258372 ,ol4 , 3656759742100547860823 ^16 , n( qlS\ ('P^r\\ 
4104021796875 ^ 150104597220703125 ~'~ 3752614930517578125 ^ t '^iM J \^ 

4^-1 _/|i52/0_316/o3| 15296 a5 6£ ,06 _ 6651852 , 1088 aS , 4502031416 /o9 

3/^ '15'- 175'- ' 2625 I- 125 336875 I- ' 4375 ~'~ 65690625 1^ 

203936 /olO _ 192671527632 All 1 , 1013665088 ,-ol2 , 1306699165544048 /ol3 _ 1577414046048 /ol4 
459375/^ 766390625 294765625 ^ 1368007265625 134118359375 

187249702135322630044 /q15 , 27983550800928 /0I6 , 197406026102200399712 _ 561057571773872512 glS 

50034865740234375 1^ ' 670591796875 >- ^ 13167069931640625 >- 4389023310546875 ^ 

9420235976547662005883368 old , 17183199427949180266496 o20 , 3314631986901973406398391696 /j2l 

154004373521044921875 ^ ' 43780507522705078125 ^ ~'~ 13090371749288818359375 
1688528041491767277609504 ,o22 _ 5344697433752690539376157188912 q23 

1328008728188720703125 5039793123476195068359375 
710002489923384145123709548576 /o24 , /nf ,Q25N fFIpl 

176392759321666827392578125 I \^ 

3/J _ 27/32 ^ 1176 ^3 _ 10233 _|_ 623418 _ 27129006 _|_ 9444396 q7 _ 82196703^8 

983643153048 /39 _ 299630551162554^10 , 169503983781937692 oil _ 1475232154782748902 q12 

9625 336875 21896875 21896875 

987637288715664984 ,jl3 _ 50792385588434618148 oli , 82664811252077369258388 olb 

1684375 9953125 1861234375 

35253092409175923838220943 0I6 , 15694806775998694225868699994 /3l7 

91200481375 ^ 1665255516875 

1993179 i277023()093KT()01i;-i2573158 0I8 , 1700015701()9775B797K7297711329632 /ol9 

68074647265625 ^ ^ 667131543203125 ''^ 

7397813830453458326408000754198846 ^20 , 7404260658028457836903051615667215212 ^21 

3335657716015625 383600637341796875 
1 172825899888916586985509963759764004972 ,o22 , 1275921695053691221896293444849460211245084 ,o23 

6981531599620703125 872691449952587890625 
1586377139730021695207597717642888690409438 ^24 , 96646235832092490263933098748637598311624744 ^25 

124670207136083984375 872691449952587890625 1^ 

221351162997126072386711604146938113185602308 q2& 

229655644724365234375 
408671731065098636806177681425625771659960813158664 q27 

48718000193603218994140625 1^ 
273870833487071174807420156540338331788208519364045012 ,o28 

3751286014907447862548828125 1^ 
5689559760891835362935463361457661566503603613628700393824 yo29 

8954319717584078047904052734375 1^ 
4208991513353764866716270905427277208225624532617948200865512 oZQ , /nf./^^lx (Vl^\ 
761117175994646634071844482421875 1^ "r ) \^ ) 



2. N ^4 
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I 909455600757890632 fj2l _ 4903884524065528 o22 _ 85272481476107112904 ,o23 , 19055775959948008 q24 

2327379429375 6630710625 48874968016875 ~'~ 5425126875 

I 821836927190253245168 o25 _ 7802995733962068568 ^26 _ 15607345431297974660104 ^27 

104732074321875 ^ 465475885875 439874712151875 

I 19608530438150678273776 ^28 , n(q'29\ 
244374840084375 t'^IM J 

_ 2 /0-2 _ 22 I 8/o2 I 494 /o4 _ 8284 ,06 , 161257/^8 _ 3197867 qW , 40478257/^12 
-'"v 9/^ 27 ' ' 1215 3645 ' 25515 ^ 127575 382725 

_ 23351747878 oU , 2016622630922 ,ol6 , f}/ a^8\ (FQc]) 
51667875 ^ ^ 1023023925 ^ ' '^{y J 
T\yr2 _ 4/0-1 _AA_32f)_56f)3, 2764 /o5 _ 16 /06 _ 30376 7 1 64 /^S _|_ 2357776 /o9 184 ol 
iWjj 'I -|- g pi "T 405 y 27!^ 1215 ~'~ 81 ~'~ 25515 135*^ 



?10 



19782236 /Oil , 76984 ol2 , 1686885896 ol3 _ 2132936 ^14 _ 5039762128^15 , 76496944^16 



54675 10935 1148175 76545 ^ 820125 ^ 688905 

'87093529676 ql7 _ 22231013972 o 18 _ 505599970587' 
113669325 ^ 51667875 ^ 442047375 



+ 



775676060228 q20 , 391764598151298668 ,o21 _ 512747827960696 o22 _ 110242538271813248108 ,o23 
465010875 ' ^ 775793143125 y 76726794375 y 48874968016875 y 

I 43888147443715136 A-?24 , /n/^/^25^ CPOo 
1611262681875 '-^KP ) ^« 

u = 3/3- 27/32 + 235/33 - 2043/3^ + 17765/3^ - 154479/3^ + 2oi4p3^7 _ 584|373 ^8 

I 4570670527 q9 _ 13248238489 10 , 7257685461779 ^11 _ 4207327698821 12 , 548778458844239 ql3 

' 45 y 15 y ' 945 y 63 y ' 945 y 

_ 3408540645268981 ql4 , 1867275863842718053 ol5 _ 51546286924640369 ol6 , 235318541138701614541 ql7 

675 ^ 42525 135 ^ 70875 

_ 61387()(i35()2S2887!)()2919 ol8 , 5284592580750019875772771 ,ol9 _ 5105851777719127773522979 ,q20 
5>15>fi25 y 5>1 049875 y 5338875 P 



212625 ^ ^ 21049875 ^ 2338875 

L53799097820463863 /o21 _ 24322483746219174453 
21049875 y 147349125 



21049875 y 147349125 ^ ' 

123726647366907086076547668505151 io23 _ 71725146330946050495871958406491 q 

86199238125 y 5746615875 P 

21829256518148380793466834521759833 ,o25 _ 569455321970865221978936851386596153 q 
201131555625 P 603394666875 P 

222828956715504343457272488503713477607 q27 _ 13181163564610119027733292785978467713 q 
27152760009375 P 184712653125 P 

16848860749534433639413131084425510515541 yo29 _ 775190057499031799755350461405515599727 ^30 
27152760009375 P 143665396875 P 



)24 

i26 

i28 



+ 0(/33i) (F2f) 



3. iV = 8 



P— q_iq3_i_46f>5_ 7552 q7 I 2302 q9 _ 555317288 ^11 , 2746333288^13 _ 1785780786208 ql5 
^ P 5P ^ 15P 525 35 1771875 P 1771875 P 227390625 P 

I 49727658484666 q\7 _ 422283776631137752 ,ol9 , 26014686098048746528 q2l 

1227909375 P 1995352734375 P 23279115234375 P 

_ 151047237869123839244528 q23 , 10528747447001620271543026 q25 _ 54722318647021454789447096 q27 

25383189111328125 P ~'~ 329981458447265625 P 318009342041015625 P 

I 141254571303908072786177902882 ,029 , n( q3l\ 
"T 151461489427294921875 P ~^ ^\P ) \r OA) 

= 1 + 3;^ + 6/32 + f /33 + 11 + 112^5 ^ ^^6 ^ 110^7 ^ 3^pS ^ 2|2^9 ^ IM^IO 

_ 12052904^11 _ 38292976 0I2 , 509694772 |0 13 , 35255323868 ql4 _ 201482297164 ql5 

84375 P 118125 P ~^ 590625 P ~^ 20671875 P 45478125 P 

_ 1185431950796 0I6 , 33950985399928 o 17 , 2298488652590486 0I8 _ 409115403009898784 ,ol9 

136434375 P 1461796875 P ~'~ 51162890625 P 3325587890625 P 

_ 781623089480829248 q20 , 2184027254831875996 q21 , 185758308159909008476 q22 

3325587890625 P 3325587890625 P 149651455078125 P 

_ 35281258231914454156624 r,23 _ 726617483139164562921284 r,24 , 2395898361931558337852 o25 

9999438134765625 109993819482421875 125707222265625 

I 556037939109294945044734 r,2& _ 4835166268399659267621472288 ^27 _ 763751109137093877154435748 o28 

15713402783203125 P 46747373280029296875 P 4006917709716796875 P 

I 711299012158799734199117996536 ,029 , 100262522489076708186582427306586 ,o30 , n( q3l\ 

1262179078560791015625 P 97187789049180908203125 P ^\P ) K^'^^J 

A/f ^ = l/9^1_4-|-56/^_608/o3_|_ 8544 /35 _ 56 ^6 _ 690104 n7 , 30632 08 , 895240544 r,9 _ 
^VJ-n 3/-* ^"'"15'^ 225*^ 875*^ 25*^ 16875-^ ~'~ 1875 ^ ~'" 5315625 ^ 

_ 2698048 olO _ 15337773008 oil , 78233152^12 , 15033061154264 ol3 _ 46357786976 14 

28125 P 20671875 P ^ 140625 P 4385390625 P 14765625 P 

_ 1502465010829864 ql5 , 11492930412128 q\& , 1111624153887072856 ql7 _ 3450172284172424 glS 
92093203125 y 664453125 y 13967469140625 y 36544921875 y 
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+ 



3734860553852716409584 ol9 , 8295415887259675624 o20 , 1982675233173138163625048 /d21 

9428041669921875 ' 16116310546875 ' 989944375341796875 

309125247582087653324456 ^22 _ 16910028827683018859169992 ^23 , 1691096970202616702325224 o24 

109993819482421875 ^ 1649907292236328125 ^ ~'~ 109993819482421875 

66840015065786580503988075712 ^25 _ 46345645264745356443888008 ^26 

1262179078560791015625 ^ 549969097412109375 ^ 

1206164175984894430511591451154792 /q27 , 57262390753050551283591565888 /o28 , /nf /329\ (F'^r) 
4373450507213140869140625 ~'~ 123743046917724609375 ' '^Kh' ) "^W 

— 2/Q-2 _ 34 4- 568 I 13784 n4 _ 69604 n6 , 241041191 _ 185964672997 glO 

gy 45 ~'~ 675'^ ^ 23625 16875 '^'^ ~'~ 15946875 ^ 2790703125 ^ 

_l_ 1392084587186 al2 _ 10666554935618557 ;Ol4 , 120272997266688239927 oW , r)(f^iS\ (F'^d) 
' 4385390625 6906990234375 ~'~ 15713402783203125 ''^\y J '^'^Z 

= 1/^-1 _zL-I-56/9_608/q3| 24652 /35 _ 224 /)6 _ 673136 f>7 , 10976 oS , 899887928 ad 

^ 15/^ 225^ ' 2625 ^ 375 16875 ^ ~'~ 5625 ^ ^ 5315625 ^ 
_ 480088^10 _ 5216665892 ^11 , 53139224 ^12 , 15464225123456 ;Ol3 _ 8888980504 14 

84375 y 6890625 ^ ~^ 2109375 ^ ~^ 4385390625 ^ 73828125 y 

_ 38922396533024104 gib , 5789142673744 /ol6 , 32467756338289108 gll _ 3641291779914748 /ol8 

2302330078125 ^ 9966796875 ^ ^ 391904296875 ^ 1279072265625 ^ 

_ 19467809813152568040512 nl9 , 85327706507419012316 o20 , 10330548206236154611036876 a21 

47140208349609375 ^ 6043616455078125 ~^ 4949721876708984375 ^ 

_ 117133256340256981616728 o22 _ 8000075417927232351906676 ,o23 , 1656430217477736394276208 o2i 
1649907292236328125 ^ 749957860107421875 ^ 4583075811767578125 ^ 

+ 0(/3'') (F3e) 

_ 2^ _ 27/5^ + ii^ /?-^ — 10179 qA 88357 /^^5 _ 3834927 ^6 _|_ 9320831 fj7 _ 404542143 ^8 

I 877S9K57841 q9 _ 38102 15382897 qV) , 4465019107354747 ,oll _ 430(i 188725191523 ,ol2 

875 P 4375 ^ 590625 ^ 65625 

I 336438262155902647 q12. _ 3785726352066159219 qlA , 9759896773518266783449 gib 

590625 1^ 765625 1^ 227390625 1^ 

_ 28239924187021013712337 glQ , 33093054294196142797506511 ^^17 _ 1612013627927391682265393 ,ol8 

10939578195 ^ 57491875 ^ 



+ 



75796875 ^ ^ 10232578125 ^ 57421875 

162079960512817503825568623029 ,ol9 _ 14211295974908998829041938103 ^^20 
665117578125 P 6718359375 ^ 



_l_ 427441298299484525828261030140363 ,o21 _ 25240543508916505408741986975703 q22 
23279115234375 ^ 158361328125 ^ 

I 152179939416162537931621662918853279261 q2Z _ 146775799365247069031550755820556182991 a24 
10999381 9482421 875 ^ 1 9991 tiqciAQSnzlRSTc; h' 



109993819482421875 ^ 12221535498046875 

11466649049267256355352323832429180274931 ,^25 _ 55297254904420195051495336139745677825041 ^26 

109993819482421875 ^ 61107677490234375 ^ 

367201599583278031923778483847484819042846311 q27 

46747373280029296875 ^ 
1062485144503393173493863633401017982717619013 q2S 

15582457760009765625 1^ 
57465089834508188027065801793670434250632397031 yo29 

97090698350830078125 P 
277399209545717425078049048151477437693717560565633 qZQ , n( fF^fl 
53993216138433837890625 P ^ '^KP ) \^ '->^ ) 
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FIGURES 



Figure 1 
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FIG. 1. The large-A^ limit of the ratio between Mq and the corresponding weak couphng 
asymptotic formula for the square, triangular and honeycomb lattices. 
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Figure 2 




Figure 3 
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FIG. 3. Asymptotic scaling test from the strong-coupling determinations of on the square 
lattice. We show curves of As/A2i, defined by Eqs.(|7§-|7|), for AT = 3, 4, 8 and for A^ = oo (exact). 
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Figure 4 
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FIG. 4. Asymptotic scaling test from the strong-coupling determinations of ^ on the triangular 
lattice. Curves of h^t/^2i for A?" = 3,4, 8 and for A?^ = oo (exact) are shown vs. ^g- 
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Figure 5 
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FIG. 5. Asymptotic scaling test from the strong-coupling determinations of ^ on the honey- 
comb lattice. Curves of Ah/A2i for iV = 3,4, 8 and for iV = cx) (exact) are shown vs. ^g. 
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TABLES 

TABLE L Summary of the large- calculations in 0{N) a models on the square, triangular 
and honeycomb lattices. All quantities appearing in this table have been defined in Sec. O. 





square 


triangular 


honeycomb 


c 


4 


6 


3 


Va 


1 


V^/2 


3v^/4 


t 






4i>s/c/3 


X 


4/c/3z 


A/cl3z 


4/C/32 


Ml 


2 


2 


2 


u 








E 


1 - l/(c/3) + \z 


1 - l/(c/3) + i2 


1 - l/(c/3) + \z 


M2 


= 2 (cosh/is - l) = 2 


M2 = |(cosh(V3Mt/2) - l) = -2 


Ml = |(cosh(V3Mh/2) - l) = 2 




Ml = 4(cosh(AJd/V2) - l) = 2 




= f(cosh{fMv)-l) =2(1 + 12) 




4^2 


4v^ 


4 



TABLE IL For various values of N we report the singularity closest to the origin on the square, 
triangular, and honeycomb lattices, as obtained by a Dlog-PA analsysis of the strong-coupling 
series of x {$x) ifli)^ the corresponding convergence radius of the strong-coupling 

expansion f3r ■ The errors we display are related to the spread of the results coming from different 
quasi-diagonal [l/m] Dlog-PA's using all available terms of the series, or a few less, while the 
difference between (3y. and should give an idea of the systematic error in the procedure. 



lattice 


A^ 


Px 




/3r 


square 


3 


0.590(1) ± i0.156(l) 


0.586(1) ± i0.157(l) 


0.61 




4 


0.557(10) ± i 0.226(4) 


0.555(5) ± i 0.225(5) 


0.60 




8 


0.467(4) ± i 0.298(3) 


0.467(3) ± i 0.293(1) 


0.55 




00 


0.321621. .(1 ±i) 


0.321621. .. (1 ±i) 


0.454841... 


triangular 


3 


0.3582(2) ± i 0.085(1) 


0.357(1) ± i 0.089(4) 


0.37 




4 


0.343(1) ± i 0.121(1) 


0.341(2) ± i 0.124(2) 


0.36 




8 


0.2901(1) ± i 0.1654(1) 


0.283(4) ± i 0.163(4) 


0.33 




00 


0.206711... ± iO. 181627... 


0.206711... ± i 0.181627... 


0.275169... 


honeycomb 


3 


±i 0.459(1) 


±i0.461(l) 


0.46 




4 


±i 0.4444(1) 


±i 0.445(2) 


0.44 




8 


±i0.4161(l) 


±i 0.4169(2) 


0.42 




00 


±i 0.362095... 


±i 0.362095... 


0.362095... 
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TABLE III. Analysis of the 14th order strong-couphng series of r = M^/Ml for = 3 on the 
square lattice, expressed in powers of E and [3. The first two lines report the values of the [l/m\ 
PA's and Dlog-PA's &i E = 1. The last two lines report the values of [l/m\ PA's and Dlog-PA's 
at /3 = 0.55 corresponding to ^ ~ 25. We show data for PA's and Dlog-PA's with Z + m > 11 and 
m > Z > 5. Asterisks mark defective PA's, i.e., PA's with spurious singularities close to the real 
axis for < 1 in the energy series case, or for (3 < 0.55 in the /3-series case. 







5/6 


6/6 


5/7 


6/7 


5/8 


7/7 


6/8 


5/9 


E = l 


PA 


* 


0.9965 


0.9967 


0.9955 


1.0126 


0.9980 


1.0007 


1.0120 




DLPA 


1.0002 


1.0011 


1.0023 


1.0005 


0.9995 








= 0.55 


PA 


0.9986 


0.9996 


1.0015 


1.0007 


1.0010 


1.0007 


1.0007 


1.0012 




DLPA 


0.9996 


1.0007 


0.9993 


1.0006 


0.9999 









TABLE IV. Analysis of the 16th order strong-coupling series of s = M"^ /Mq for iV = 3 on 
the square lattice. The first two lines report the values of the [l/m] PA's and Dlog-PA's at E = 1. 
The last two lines report the values of [Z/m] PA's and Dlog-PA's at /? = 0.55. We show data for 
PA's and Dlog-PA's with Z -|- m > 13 and m > Z > 5. Asterisks mark defective PA's. 







6/7 


5/8 


7/7 


6/8 


5/9 


7/8 


6/9 


5/10 


8/8 


7/9 


6/10 5/11 


E = l 


PA 


0.9947 


0.9938 


0.9941 


0.9942 


* 


0.9944 


1.0020 




0.9961 


1.0028 


1.0028 * 




DLPA 


0.9941 


0.9971 


* 


0.9992 


0.9978 


0.9951 


0.9973 


0.9984 








= 0.55 


PA 


* 


0.9971 


0.9972 


0.9974 




0.9976 


0.9988 


0.9998 


0.9980 


1.0023 


0.9996 0.9995 




DLl'A 


0.9971 


0.9980 


0.9971 


0.9992 


0.99,S5 


0.9977 


0.99.S2 


0,997() 









TABLE V. Analysis of the 20st order strong-coupling series of E-^u{E) and p~^u{/3), where 
u = m^/ix'^'i): for TV" = 3 on the square lattice. The first two lines report the values of -u as 
obtained from the [Z/m] PA's and Dlog-PA's &i E = 1. The last two lines report the values of u 
from [Z/m] PA's and Dlog-PA's at /? = 0.55. The analysis detected a pole at Eq = -0.086418 in 
the energy series, and at /3o = —0.085545 in the j3 series, corresponding to Mq = —16.000. We 
show data for PA's and Dlog-PA's with Z + m > 16 and m > Z > 8. Asterisks mark defective 
approximants. 







8/8 


8/9 


9/9 


9/10 


8/11 


10/10 


9/11 


8/12 


E = 1 


PA 


0.2491 


0.2502 


0.2495 


0.2488 


0.2491 


0.2495 


0.2504 


0.2496 




DLPA 


0.2497 


0.2524 


0.2510 


0.2486 


0.2492 








= 0.55 


PA 


0.2493 


0.2496 


0.2488 


0.2496 


0.2500 


* 


0.2495 


0.2503 




DLPA 


0.2493 


0.2495 


0.2491 


0.2498 


0.2500 
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TABLE VI. In this Table we summarize our strong-coupling results for N = 3, giving the 
estimates of r* , s* and u* from the PA and Dlog-PA analysis of both the energy and /3-series of r, 
s, and u. For all lattices considered the values of (3 where the /3-approximants have been evaluated 
correspond to a correlation length ^ ^ 20. 













s 


u 


square 


E = 


■ 1 


PA 


1.004(8) 


1.000(5) 


0.2498(6) 








DLPA 


1.0000(12) 


0.997(2) 


0.249(2) 




P = 


0.55 


PA 


1.0009(6) 


1.000(2) 


0.2499(6) 








DLPA 


1.0002(6) 


0.9978(8) 


0.2499(5) 


triangular 


E = 


: 1 


PA 




1.000(4) 


0.2497(15) 








DLPA 




0.997(3) 


0.248(2) 




P = 


0.33 


PA 




0.9985(13) 


0.2504(3) 








DLPA 




0.9980(9) 


0.2499(4) 


honeycomb 


E = 


: 1 


PA 


1.01(4) 


0.999(4) 


0.250(2) 








DLPA 


0.991(13) 


0.999(3) 


0.247(2) 




P = 


0.85 


PA 


1.001(2) 


0.9987(5) 


0.2490(3) 








DLPA 


1.0009(8) 


0.9987(5) 


0.2491(3) 



TABLE VII. The strong-coupling estimates of are compared with the available Monte Carlo 
results on the square lattice for various values of N. The strong-coupling estimates of come from 
the plain series of C^, and from ([9/10], [8/11], [9/9], [8/10]) Dlog-PA's of The iV = 3,4,8 

M.C. data are taken respectively from Refs. |p8| , p9| , [p4| . The asterisk indicates that the number 



concerns ^cxp, and not (^g- 


N (3 


plain series 


Dlog-PA's 


M.C. 


3 1.4/3 


6.567 


6.869(1) 


6.90(1)* 


0.5 


9.939 


11.036(4) 


11.05(1)* 


1.6/3 


15.429 


18.90(2) 


19.00(2)* 


1.7/3 


24.300 


33.9(1) 


34.44(6)* 


0.6 


38.459 


61.0(4) 


64.7(3)* 


4 0.45 


4.665 


4.672(1) 


4.67(1) 


0.5 


7.845 


7.87(1) 


7.83(1) 


0.55 


13.879 


13.88(5) 


13.99(3) 


0.575 


18.701 


18.6(2) 


18.91(5) 


0.6 


25.329 


24.8(4) 


25.5(2) 


8 0.5 


5.432 


5.459(1) 


5.461(5)* 


0.525 


6.584 


6.651(1) 




0.55 


7.981 


8.139(2) 




0.575 


9.659 


10.01(1) 


9.884(13)* 
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TABLE VIII. Analysis of the 11th order strong-coupling series of s = M"^ /Mq for = 3 on 
the triangular lattice. The first two lines report the values of the [l/m] PA's and Dlog-PA's at 
E = 1. The last two lines report the values of \l/m\ PA's and Dlog-PA's at /? = 0.33 corresponding 
to ^ ~ 22. Asterisks mark defective PA's. 







4/4 


4/5 


5/5 


4/6 


5/6 


4/7 


E=l 


PA 


0.9993 


0.9972 


1.0005 


0.9927 


0.9954 


1.0039 




DLPA 


0.9963 


1.0014 


0.9948 


* 






= 0.33 


PA 


0.9993 


0.9989 


1.0005 


0.9969 


0.9974 


0.9995 




DLPA 


0.9987 


* 


0.9972 


0.9975 







TABLE IX. Analysis of the 14st order strong-coupling series of E-^u{E) and P~'^u{f3), where 
u = m^/ {xmi), for = 3 on the triangular lattice. The first two lines report the values of u as 
obtained from the \l/m\ PA's and Dlog-PA's at E = 1. The last two lines report the values of u 
from [l/m] PA's and Dlog-PA's at /? = 0.33 corresponding to ~ 22. A pole has been detected at 
Eq = —0.050655, corresponding to Mq = —16.000. Asterisks mark defective PA's. 







5/6 


6/6 


5/7 


6/7 


5/8 


7/7 


6/8 


5/9 


E = l 


PA 


0.2442 


0.2502 


0.2533 


0.2483 


0.2492 


0.2494 


0.2500 


0.2497 




DLPA 


0.2433 


0.2521 


0.2502 


0.2477 


0.2491 








/3 = 0.33 


PA 


0.2521 


0.2500 


0.2500 


0.2502 


0.2502 


0.2504 


0.2505 


0.2502 




DLPA 


0.2496 


0.2494 


0.2494 


0.2496 


0.2502 









TABLE X. Analysis of the 19th order strong-coupling series of r//3, where r = M^/M^, for 
= 3 on the honeycomb lattice, expressed in powers of E and (3. The first two lines report the 
values of the [l/m] PA's and Dlog-PA's at E = 1. The last two lines report the values of [l/m] 
PA's and Dlog-PA's at /3 = 0.85 corresponding to ^ ~ 22. Asterisks mark defective PA's, i.e., PA's 
with spurious singularities close to the real axis for < 1 in the energy series case, or for (3 < 0.85 
in the /^-series case. 







8/8 


7/9 


8/9 


7/10 


9/9 


8/10 


7/11 


9/10 


8/11 


7/12 


E 


= 1 PA 


1.070 


0.963 


1.006 


1.035 


* 


0.981 


1.035 


0.980 


0.980 


1.080 




DLPA 


1.007 


0.996 


0.977 


1.010 


0.993 


0.989 


* 










= 0.85 PA 


1.0039 


* 


1.0020 


1.0031 


* 


0.9988 


* 


0.9996 


0.9997 


1.0024 




DLPA 


1.0006 


0.9992 


1.0019 


1.0009 


1.0009 


1.0014 


1.0005 










TABLE XL 


Analysis of the 25th order strong 


5-couplin 


g series 


of S = 




for N 


= 3 on 



the honeycomb lattice. The first two lines report the values of the [l/m] PA's and Dlog-PA's at 
E = 1. The last two lines report the values of [l/m] PA's and Dlog-PA's at /3 = 0.85 corresponding 
to ^ ~ 22. Asterisks mark defective PA's. 







11/11 


10/12 


11/12 


10/13 


12/12 


11/13 


10/14 


12/13 


11/14 


10/15 


E=l 


PA 


0.9956 


1.0001 


1.0052 


0.9983 


* 


* 


* 


* 


* 


* 




DLPA 




0.9964 


0.9994 


0.9963 


1.0023 


* 


0.9972 








= 0.85 


PA 


0.9978 


0.9984 


0.9983 


0.9982 


* 


0.9983 


0.9991 


0.9982 


0.9989 


0.9989 




DLPA 


0.9989 


0.9982 


* 


0.9982 


* 


0.9992 


0.9983 
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TABLE XII. Analysis of the 29st order strong-coupling series of E-^u{E) and where 
u = m2/(xm4), for = 3 on the honeycomb lattice. The first two lines report the values of u as 
obtained from the [l/m] PA's and Dlog-PA's at E = 1. The last two lines report the values of u 
from [l/m] PA's and Dlog-PA's at /3 = 0.85 corresponding to ^ ~ 22. A pole has been detected at 

Eq = -0.11101. corn-spoudiug lo = -IG.OOO. AsU-risks mark defective PA's. 







13/13 


13/14 


12/15 


14/14 


13/15 


12/16 


14/15 


13/16 


12/17 


E = 1 




0.2511 


0.2181 


0.2198 






* 


0.2182 


0.2195 


0.2526 




DLl'A 


().2i()5 


0.211)9 


0.218."') 


0.2107 


0.211)7 










P = 0.85 


PA 


0.2490 


0.2490 


0.2490 


0.2489 


0.2482 


0.2488 


0.2489 


0.2490 


0.2490 




DLPA 


0.2493 


0.2494 


0.2488 


* 


* 


* 









TABLE XIII. Singularities in the complex /3-plane for the triangular and honeycomb lattice 
at A'' = oo with positive real and imaginary part for the lowest values of (d, c). The singularity on 
the the real axis for the honeycomb lattice (/? = 0.627168) does not appear on the principal sheet 
of w{P) as the corresponding tw-value is 0.962998i. 



{d,c) 


triangular 


honeycomb 


(1,0) 




0.3620955333 i 


(1,±2) 


0.206711 + 0.181628 i 
0.685669 + 0.749077 i 


0.482696 + 0.628020 i 
0.449772 + 0.583632 i 
0.627168 


(3, ±4) 


0.240692 + 0.486530 i 
0.564118 + 0.203430 i 
1.469137 + 2.118380 i 


0.566020 + 1.476842 i 
0.946032 + 1.663513 i 
1.237526 + 0.266631 i 


(5, ±6) 


0.260362 + 0.780732 i 
0.920774 + 0.210433 i 
2.244964 + 3.482094 i 


0.627495 + 2.353352 i 
1.413547 + 2.691806 i 
1.836347 + 0.535158 i 


(5, ±8) 


0.662428 + 0.858005 i 
0.980284 + 0.576976 i 
2.847075 + 3.621622 i 


1.502252 + 2.655737 i 
1.907210 + 2.926524 i 

2.022490 + 1.997953 i 
2.500324 + 0.265750 i 


(7, ±8) 


0.274172 + 1.072281 i 
1.763190 + 0.911108 i 
3.019312 + 4.844833 i 


0.669817 + 3.225087 i 
1.881646 + 3.718975 i 
2.432118 + 0.804519 i 
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TABLE XIV. For N = oo we report the zeroes closest to the origin as obtained by an analysis 
of the strong-couphng series of x by Dlog-PA's and lA's. We consider the series at 15th and 
30th order on the triangular lattice, and 30th and 60th order on the honeycomb lattice. 7^ is the 
exponent corresponding to the singularity Ps in the lA analysis [^0| (its exact value is 7^ = — 1 /2) . 
The Dlog-PA analysis does not provide stable estimates of 7^ . The values we quote are average and 
maximum spread of the Dlog-PA's [m/n] with 5 < m < 9 (series with 15 terms), 12 < m < 17 (30 
terms) and 27 < m < 32 (60 terms); for the lA we use in all cases 6 quasi-diagonal approximants. 



lattice 


n 


Ps (Dlog-PA) 


Ps (lA) 


7s (lA) 


triangular 


15 
30 


0.214(4) ±i0.1838(2) 
0.2084(2) ±iO. 1821(2) 


0.206(1) ±i 0.1811(6) 
0.206712(1) ±i 0.181628(1) 


0.54(7) =F« 0.07(6) 
-0.4997(4) =F« 0.0001(2) 


honeycomb 


30 
60 


±i 0.3648(8) 
±i 0.36270(4) 


±i 0.36211(3) 
±i 0.3620955327(5) 


-0.50(2) =F« 0.001(6) 
-0.500001(1) 0.000001(1) 
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